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O ; ABSTRACT 

i We examine the growth of the stellar content of galaxies from z = 3 ^ in cosmologi- 

• ' cal hydrodynamic simulations incorporating parameterised galactic outflows. Without 

! outflows, galaxies overproduce stellar masses {A'h) and star formation rates (SFRs) 

'p^' compared to observations. Winds introduce a three-tier form for the galaxy stellar 

Q I mass and star formation rate functions, where the middle tier depends on differential 

, (i.e. mass-dependent) recycling of ejected wind material back into galaxies. A tight 

C/5 ' M*— SFR relation is a generic outcome of all these simulations, and its evolution is 

well-described as being powered by cold accretion, although current observations at 

z > 2 suggest that star formation in small early galaxies must be highly suppressed. 

' Roughly one-third of z = galaxies at masses below M* are satellites, and star forma- 

^ , tion in satellites is not much burstier than in centrals. All models fail to suppress star 

OO ' formation and stellar mass growth in massive galaxies at z < 2, indicating the need 

CNJ ' for an external quenching mechanism such as black hole feedback. All models also fail 

, to produce dwarfs as young and rapidly star-forming as observed. An outflow model 

CO ■ following scalings expected for momentum-driven winds broadly matches observed 

\ galaxy evolution around M* from z = — 3, which is a significant success since these 

, galaxies dominate cosmic star formation, but the failures at higher and lower masses 

T-H ' highlight the challenges still faced by this class of models. We argue that central star- 

I forming galaxies are well-described as living in a slowly-evolving equilibrium between 

■ inflows from gravity and recycled winds, star formation, and strong and ubiquitous 

• i-H ' outflows that regulate how much inflow forms into stars. Star-forming galaxy evolution 

, is thus primarily governed by the continual cycling of baryons between galaxies and 

H ' intergalactic gas. 
■ 

Key words: galaxies: formation, assembly, low redshift; intergalactic medium; cos- 
mology: theory; methods: numerical 



1 INTRODUCTION 

Galaxies are the primary tracers by which astronomers map 
and measure the Universe, so understanding how galaxies 
form and evolve is a longstanding central goal of astronomy. 
The physics of galaxy formation involves a diverse set of 
processes spanning an enormous dynamic range, from black 
holes and stellar physics on sub-parsec scales to cosmology 
and large-scale structure on gigaparsec scales. As a conse- 
quence, despite rapidly advancing observations of galaxies 
across cosmic time, a comprehensive theory for the forma- 
tion and evolution of gala^xies that can explain all their ob- 
served properties remains elusive. 

As galaxies are essentially collections of stars, under- 
standing the growth of the stellar component of galaxies 



is of fundamental importance. Recent observations have 
probed the star formation rates and stellar masses of galax- 
ies out to redshift z ~ 7. Key trends have emerged. For in- 
stance, fairly massive galaxies are in place at early epochs, 
and galaxies exhibit "downsizing" in the sense that more 
massive galaxies form their stars earlier and have lower 
birthrates today. While some claim that these trends are in 
conflict with expectations from hierarchical structure forma- 
tion, they are actually not. Since the most massive galaxies 
form within the largest mass density peaks, they collapse 
first a nd start forming stars earlier than lower-mass sys- 
tems I Dave et al.l l2006l : iNeistein. van den Bosch, fc DekeTI 
'2006^). The collapse is driven by gravitational instability and 
strong radiative cooling that (left unchecked) grow galax- 
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ies very rapidly at early epochs, so that models without 
strong fee dback yield galaxies that are too large at high red- 
shifts (e.g. IPave. Finlator. fc Oppenheimer 1120061 ). It is im- 
portant to note that no confirmed high-z galaxy yet observed 
is too large to be produced in a hierarchical model, and in 
fact strong feedback appears to be necessary to suppress 
gala:xy formation from the earliest epochs. The evolution of 
galaxy stellar masses and star formation rates therefore pro- 
vide valuable insights into the processes, particularly related 
to feedback, that drive galaxy formation from the highest 
redshift until today. 

The canonical scenario for galaxy formation posits 
that galaxies assemble first as rotationally-supported disks 
cooling out of virial-temperature ga s within a dark mat- 
ter halo (e.g. iRees fc Qstrikeil 1 19771 : IWhite fc Reei 1 19781 : 
IWhite fc FrenkI 1991 ) , and then such objects merge to form 
large r and earlier-type systems (e.g. IMih os & Hernguistl 
I1996D . Naively, one expects then that star- forming disks 
would be surrounded by hot gaseous halos (| Benson et al.l 
I2OOOI ). and that early low-mass objects would tend to be 
disky while more massive objects would be ellipticals. Re- 
cent observations, however, do not clearly support this 
paradigm. For instance, star-forming disks and hot (X-ray 
emitting) gaseous halo s are rarely if ever found toget her (e.g. 
iRasmussen et al"]|2009l . though see lCrain et al.|[2010l and ref- 
erences therein) . Observed dynamics of 2; ~ 1 — 2 galaxies in- 
dicate that larger objects te nd to be disky, while smal l er one s 
are dispersion-dominated (jForster Schreiber et al.l l2009f ). 
Furthermore, galaxies appear to be consuming their 
gas rapidl y enough that fue l must be constantl y re- 
plenished ((Tacconi et all I2OI0I : [Papovich et al] I2OI0I '). but 
merger-induced starb ursts do not appears to be driving 
global star formation |Brinchmann et al ]|2004l : lNoeske et al.l 
l2007l : |jogee et aLlbOOgf T Hence there is some doubt that the 
canonical halo/merger-centric view provides the appropriate 
framework for galaxy formation, despite the wide-ranging 
successes tha t models based on this have enjoyed (e.g. 
lBenson|[2010l ). 

Improvements in cosmological hydrodynamic simula- 
tions over the past decade have made them a competitive ap- 
proach towards understanding the physics of galaxy forma- 
tion in a hierarchical structure formation context. The com- 
plex three-dimensional geometry of large-scale structure is 
self-consistently evolved together with the baryonic compo- 
nents that make up observable galaxies and the co-evolving 
intergalactic medium (IGM). These simulations have long 
suggested that accretion onto galaxies is powered by cold, 
filamentary streams that connect to the cosm ic web on 
larger scales (|Keres et al.ll2005l : [Pekel et al.ll2009l '). and that 
most accreted gas f ueling star formation never heats to the 
virial temperature (|Katz fc Gunnlll99ll :[ Dekel fc BirnboimI 
I2OO6I ). Gas replenishment from the IGM is a natural out- 
come of this model, and generally occu rs smoothly rather 
than being driven by major mergers (e.g. Murali et al.ll200"3 : 
iKeres et al] l2005l : Ivan de Voort et al.l bOlOD . These simula- 
tions provide many insights into how galaxies form and grow, 
some of which are in common with the traditional scenario, 
while some differ. 

Recently, simulations incorporating strong and ubiq- 
uitous outflows have had success matching a wide range 
of data on galaxies (e.g. iDave. Finlator. fc Oppenheimer I 



the IGM (e.g. lOppenheimer fc Davi I2OO6I . I2OO8I . l2009l : 
lOppenheimer. Dave, fc FinlatoJ \200% . The outflow mod- 
els in these simulations were originally tuned to enrich 
the IGM in a ccord with observations o f 2: ~ 2 — 4 
C IV absorbers jOppenheimer fc DavdlioO^ ). The favoured 
scalings turn out to be those expected for mom e ntum- 
conserved winds (Murray. Quataert. fc ThompsonI l2005l : 



'Zhang fc ThompsonI l20ld ). namely that the wmd speed 
scales with circular (or escape) velocity, while the amount 
of mass ejected per unit star formation (i.e. the mass 
loading factor) scales inversely with it. This provides suf- 
ficiently early IGM enrichment by expelling the major- 
ity of metals formed in small high-redshift galaxies, while 
having moderate wind speeds that does n ot to over- 
heat the IGM lOppenheimer fc Pavel l2006t ). In subse- 
quent work we showed that an outflow model with sim- 
ilar scalings consistently does as well or better than any 
other outflow model that we have tried at matching a 
wide range of IGM properties. We note that the ac- 
tual implementation of this wind model in our simula- 
tions has evolved owing to advances in modeling techniques 
and changes in cosmol ogy (fOpp cnhcimcr fc Dave 200^; 
lOppenheimer et aL 20101'). Successes include very earl y {z > 



, ery 1 . , 

5) enrichment l Oppenheimer. Dave, fc Finlatoill2009l ). Ovi 
absorbers at low-z (with the introd uction of a physically- 
moti vated IGM turbulence model lOppenheimer fc Pavel 
i2009i ) , the metallicities in intragroup and intracluster 
gas (|Dave. Oppenheimer. fc SivanandamI I2OO8I ). and the 
overall cosmic metal distr ibution (including the so-c alled 
"missing metals problem" : iDave fc Oppenheimeill2007l 'l. 

Examining galaxy properties, we again found that an 
outflow model following similar scalings was the most suc- 
cessful amongst the ones we tried. Successe s included match- 
ing the mass-metallicity relation at 2: = 2 l|Finlator fc Davi 
I2OO8), t he shape of the pres ent-day stellar mass func- 
tion ^Oppenheimer et al. 2OI0I) below M* , damped Lya 
absorber kinematics ( Hong et al.l |2010|). and observations 
of very high-g (z > 4) galaxies ( Finlator et al.l 2006 ; 



I2OO6I : iFinlator fc DavelboOSi : lOppenheimer et al.ll2010l ) and 



^Dave. Finlator. fc Oppenheimer II2OO6I : Bouwens et al.l 2007 ; 
Finlator et al. 201(|). This is not to say that this par- 
ticular outfiow model matches all the data; in fact, it 
does not, it simply d oes better than other s attempted, in- 
cluding the oft-used ISpringel fc Hernguistl (|2003bl ) model. 
Moreover, the implementation remains highly heuristic, and 
m ay be impacted by numer ical effects (e.g. see discussion 
in lOppenheimer et al.ll2010l ). Nevertheless, the consistent 
(comparative) success of an outfiow model with momentum- 
driven wind scalings suggests that it captures some essential 
ingredients required for a successful model of galactic out- 
flows, and motivates a more detailed investigation of the 
implications of such outflows. 

Two key predictions from models incorporating out- 
flows with such scalings are that (i) the total amount of 
material eject ed from galaxies signiflcantly exceeds that go- 
ing into stars (|Oppenheimer fc Dave|[200a hereafter OD08); 
and (ii) the ejected material often re-accretes onto galax- 
ies lOppenheimer et al.ll2010l '). The strong and ubiquitous 
ejection and re-accretion of material suggest that galaxies 
are best viewed as evolving within an ecosystem of surround- 
ing intergalactic gas with which they continually exchange 
mass, energy, and metals. 

While the above results are encouraging, these simu- 
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lations have yet to be comprehensively tested against ob- 
servations of galaxy evolution from the present epoch back 
into the early universe. Such a comparison would provide the 
most stringent test to date of how well these simulations can 
reproduce the observed Universe, and would better highlight 
the current successes and failures of galaxy evolution simu- 
lations. This series of papers aims to understand how basic 
galaxy statistics and scaling relations between their stellar, 
gas, and metal content constrain the physics of galaxy for- 
mation. In Paper I (this paper) we focus on the statistics of 
and scaling relations between stellar mass and star formation 
rate. In Paper II (Dave, Finlator, & Oppenheimer, in prep.) 
we will additionally investigate the gas and metal content of 
galaxies. The primary goal of these works is to understand 
and characterise the key physical processes that drive the 
observable properties of galaxies across cosmic time. 

In these papers, we will demonstrate that for a par- 
ticular choice of outflow model, namely outflows following 
scalings as expected for momentum-driven winds, the sim- 
ulations can reproduce key observations of scaling relations 
between stars, gas, and metals for large star-forming galax- 
ies that dominate cosmic star formation. This is a first for 
cosmological hydrodynamic simulations. However, all simu- 
lations we examine here fail to match selected galaxy prop- 
erties at high and low masses, indicating that additional 
or different physics is required even for our most success- 
ful model. Using insights gained from the simulations, we 
show that relations between basic galaxy constituents can 
be broadly understood through a simple scenario in which 
inflows are balanced by outflows and star formation. This 
contrasts with the canonical halo/merger-centric model of 
galaxy formation; in our scenario, the dark matter halo's 
virial radius and galaxy mergers play a sub-dominant role in 
the evolution of star-forming systems. We argue that galac- 
tic outflows are the key moderator of the stellar content, star 
formation rate, and gas and metal evolution in galaxies at 
all epochs. In turn, observations of these properties provide 
valuable constraints on the cosmic ecosystem within which 
galaxies form and grow. 

This paper is organised as follows. In ^ we describe 
our hydrodynamic simulations including our galactic outflow 
model. In 33] we investigate the galaxy stellar mass and star 
formation rate functions, and explain the three-tier behavior 
generically seen in wind simulations by invoking differential 
wind recycling. In 33]we look at the scaling relation between 
SFR and M« and its evolution to high-z, emphasizing its 
connection with the matter accretion rate into halos. In ^5] 
we further break down these trends in terms of centrals vs. 
satellite galaxies, and raise some issues regarding observed 
dwarf galaxies that none of the models reproduce. Finally, 
we summarise and discuss the broader implications of our 
work in [JS] 



2 SIMULATIONS 
2.1 Code 

Our simulations are run with an extended version of the 
Gadget-2 N-body -I- Smo othed Particle Hydrodynamic 
(SPH) code l|Springel|[2005l . OD08). We assume a ACDM 
cosmology using the c osmological paramete rs based on re- 
cent WMAP results jHinshaw et all 120091 '): = 0.28, 



Qa = 0.72, h = Ho/{WO kms"^ Mpc"^) = 0.7, a primordial 
power spectrum index n = 0.95, an amplitude of the mass 
fluctuations scaled to erg — 0.82, and fib — 0.046. We refer to 
these runs as the r-series, where our general naming conven- 
tion for a simulation run is T[boxsize]n[particles/side]['wind 
model]. Our primary runs use the boxsize of 48/i~^Mpc on 
a side with 384^ dark matter and 384^ gas particles. To 
expand our dynamic range we also run two additional sets 
of simulations that are identical except that they have box 
sizes of 24/i~^Mpc and 48/i~^Mpc, and have only 256^ par- 
ticles of each species. The particle masses and gravitational 
softening lengths are listed in Table [T] 

An overview of the Gadget-2 code can be found 
in §2 of K09a. Additions to the public version of the 
code include cooling pr ocesses usin g the p rimordial abun- 
dances as described in 'Katz et al.' (l99&) and metal-line 
cooling as described in £)ppcnhcimor & Dave (2006). Star 
formation is modeled using a subgrid recipe introduced in 
ISpringel fc Hernguistl l|2003h where a gas particle above a 
critical density is modeled as a fraction of cold clouds embed- 
ded in a warm ionised medium following iMcKee fc Ostrikei 
ilQTlh . Star formation (SF) foUows a Schmidt law (jSchmidt 
I1959I ) with the SF timescale scaled to match the z = Ken- 
nicutt law (lKennicuttlll998al '). The dens it y thr eshold for SF 
is riH = 0.13 cm~ . We use a lChabrierl (|2003l ) initial mass 
function (IMF) throughout. The fraction of mass of the IMF 
going into massive stars, defined here as ^ 10 Mq and as- 
sumed to result in Type II supernovae (SNe), is /sn = 0.18; 
for reference, a Salpeter IMF has /sn = 0.10. 

Our simulations directly account for metal enrichment 
from sources including Type II supernovae (SNe), Type la 
SNe, and asymptotic giant branch (AGB) stars. Gas par- 
ticles eligible for SF undergo self-enrichment from Type 
II SNe using the instantaneous recycling approximation, 
where mass, energy, and metallicity are assumed to im- 
mediately return to the ISM. Type II SN metal enrich- 
me nt uses the metallic it y-dep endent yields calculated from 
the IChieffi fc Limongil (120041 ) SNe models. Our prescrip- 
tions for feedback, described below, also assume the in- 
stantaneous input of energy from O and B stars (i.e. stars 
th at are part of /sn). W e inpu t the Type la SNe rates 
of IScannapieco fc Bildster] (|2005l ), where an instantaneous 
component is tied to the SFR, and a delayed component 
is proportional to the stellar mass as described in OD08. 
Each Type la SN adds 10 ^^ ergs of thermal energy and the 
calculated metal yields of [Thielemann et al.l (|l986l ) to sur- 
rounding gas; the mass returned is small compared to that 
from Type II SNe. AGB stellar enrichment occurs on delayed 
timescales from 30 Myrs to 14 Gyrs, using a star particle's 
age and metallicity to calculate mass loss rates and metallic- 
ity yields as described in OD08. AGB mass and metal loss 
is returned to the three nearest surrounding gas particles. 
OD08 showed that the largest effect of AGB stars is to re- 
plenish the ISM, because a star particle can lose over 30% 
of its mass over tu, nearly double the 1 — /sn = 18% that is 
recycled instantaneously via Type II SNe. 

Galactic outflows are implemented using a Monte Carlo 
approach analogous to star formation. Outflows are directly 
tied to the SFR, using the relation Mwind = V^sf, where 
is the mass loading factor. The probabilities for a gas parti- 
cle to spawn a star particle is calculated as described above, 
and the probability to be launched in a wind is r; times that. 
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If the particle is selected to be launched, it is given an ad- 
ditional velocity of in the direction of vxa, where v and 
a are the particle's instantaneous velocity and acceleration, 
respectively; this would create a bipolar outflow for a thin 
disk, but in practice the outflow has a much larger opening 
angle. This subgrid model circumvents our inability to re- 
solve the detailed physics of the ISM that generates winds, 
at the cost of introducing additional parameters. Choices of 
the parameters i] and define the "wind model" . 

Once a gas particle is launched, its hydrodynamic 
(not gravitational) forces are turned off until either 1.95 x 
10^'^/(vw( kms~^)) years have passed or, more often, the gas 
particle has reached 10% of the SF critical density. This is 
intended to simulate the formation of hot chimneys extend- 
ing out of a disk galaxy providing a low resistance avenue 
for outflows to escape into the IGM, which the spherically- 
averaged SPH algorithm at our current resolution is inca- 
pable of resolving. However, explicit decoupling neglects en- 
ergy losses owing to the initial creation of such channels as 
well as any drag along their interfaces, and as such cannot 
drive turbulence in the ISM. We note that decoupling has a 
significant impact on wind propagat ion; simulations without 
such decoupling explored e.g. by iDalla Vecchia fc Schavd 
l|2008t l show that decoupling greatly enhances the ability of 
material to escape a galaxy's ISM. Hence decoupling must 
be regarded as an integral part of the wind models described 
in this work, having a significant impact on the results, par- 
ticularly those associated with recycling of wind material 
(which would be much more prevalent if the hydrodynamic 
forces were never turned offl. lSpringel fc Hernguisl l|2003bl ') 
demonstrated that this decoupling achieves resolution con- 
vergence in the cosmic star formation history. We will later 
explicitly demonstrate resolution convergence for many of 
the galaxy properties of interest here. 

2.2 Wind Models 

For this paper we run four wind models, described as follows: 
No winds (nw): We turn off our galactic winds. Note 
that, as with all simulations, energy is imparted thermally 
to ISM SPH particles using the lSpringel fc Hernguisl (|2003l ) 
subgrid two-phase recipe where all SN energy is instanta- 
neously coupled to the ISM. Owing to the explicit coupling 
of the hot and cold ISM phases, such energy cannot drive a 
wind. 

Constant winds (cw): Feedback energy is added ki- 
netically to gas particles at a constant rate relative to the 
SFR, with 77 = 2 and a constant velocity w„ = 680 kms~^. 
This is intended to mimic the constant-wind model of 
ISpringel fc Hernguistl (|2003l ). but with some notable differ- 
ences. Here r] is defined relative to the SFR of the entire 
Chabrier IMF, and not just the long-lived stars in a Salpeter 
IMF as it is defined by^pringcl & Hornquist (2003). Second, 
Ww = 680 kms~^ is used rather than 484 kms~^, because 
there is more SN energy per mass formed in a Chabrier IMF. 
The kinetic energy imparted to wind particles per unit SF 
is 9.25 X 10"*** erg/ M©, which is 95% of the SN energy in this 
IMF assuming aU stars ^ 10 Mq add 10^^ ergs/SN. The 
hydrodynamic decoupling time in this model has a maxi- 
mum of 2.9 X 10'' years. We note that this model is similar 
to the "REFERENCE" outfiow model of the Over Whelm- 
ingly Large Simulations (OWLS ISchave et al.ll2010l ) , except 



that we include hydrodynamic decoupling while they do not; 
in fact, they have also considered models likw nw and vzw 
(below) except without decoupling. 

Slow winds (sw): We run an alternative constant- 
wind model with the only difference being the outflow ve- 
locities are half as high as the cw winds (i.e. i;„ = 340 kms~^ 
and rj — 2). Therefore only roughly a quarter of the SN en- 
ergy couples kinetically into these outflows. 

Momentum-conserved winds (vzw): 

This model uses the a na lytic derivations by 
[Murray. Quataert. fc ThompsonI (|2005l ) based on the 
galaxy veloci ty dispersion (a). I ts imp lementation is 
explained in lOppenheimer fc Pavel (|2006l) and updated 
in OD08 with the addition of an in s itu group flnder to 
calcul ate a using galaxy mass following IMo. Mao, fc White! 
(1998!). The wind parameters vary between galaxies using 
the following relations: 



V ^ —, 

a 



(1) 
(2) 



where /l is the luminosity factor in units of the galactic Ed- 
dington luminosity (i.e. the critical luminosity necessary to 
expel gas from the galaxy potential), and ctq is the normal- 
ization of the mass loading factor. Here we randomly select 
/_L — [1.05, 2] for each SPH particle, following observations 
bv lRupke. Veilleux fc Sanders! (|2005l ). and include a metal- 
licity dependence for /l owing to more UV photons output 
by lower metallicity stellar populations 



h = /l,© X 10 



-0.0029*{log Z+9y- ^+0.4,17694 



(3) 



We further add an extra kick to get out of the potential of 
the galaxy in order to simulate continuous pumping of winds 
until it is far from the galaxy disk (as argued by MQT05). 
These parameters and /l are well-constrained by obser- 
vations, leaving the only unconstrained free parameter as ctq. 
We set (Jo = 150 kms~^, which is the same efficiency factor 
of the mass loading used in OD08. The energy budget for 
winds scales as cx rjv'^ oc a, and exceeds the supernova en- 
ergy budget only at quite large halo masses of > 10^* Mq (at 
z = 0). Since we do not include any type of f eedback that 
quenches star formation in massive galaxies (|Gabor et al.l 
2010), our simulations do not accurately model these sys- 
tems in any case. In general, the wind energetics in this 
model are quite modest, and despite the large mass outfiow 
rates and speeds around the escape velocity, are generally 
well below the supernova energy budget (globally, about half 
in the r48n384vzw run). This does not even account for ra- 
diative energy input from massive stars, which can exceed 
the energy output from supernovae. 

These wind models allow us to illustrate how out- 
fiows impact the scaling relations of galaxies. An impor- 
tant difference comes from the wind speed scaling: In the 
momentum-conserved outflows case, the wind speed is typ- 
ically around the galaxy escape velocity, and scales with 
it. Hence in this model, the potential well of the galaxy 
does not play a governing role in wind dynamics as a func- 
tion of mass. This stands in contrast to the canonical model 
for outflows (e.g. iDekel fc Silklll986! ) in which outflows es- 
cape from small galaxies owing to their small potential 
wells while being recaptured in large galaxies. If gravity 
was the only relevant force, this model should show ex- 
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Table 1. Simulation runs 



Name 




e'' 


mspH 


^dark 


gal 


( km s ^ ) 




wiii*_i/ oiN 


Pwind/P» 




Wind 


r48n384nw 


48 


2.5 


3.6 X 10^ 


1.8 X 10* 


-lOlo 














0.220 





r48n384vzw 


48 


2.5 


3.6 X 10^ 


1.8 X 10* 


1.1 X 10^ 


OC (T 


OC cr~^ 


0.49» 


5.2 


0.097 


0.345 


r48n384cw 


48 


2.5 


3.6 X 10^ 


1.8 X 10* 


1.1 X 10^ 


680 


2 


0.95 


8.7 


0.045 


0.174 


r48n384sw 


48 


2.5 


3.6 X 10^ 


1.8 X 10* 


1.1 X 10^ 


340 


2 


0.24 


4.4 


0.124 


0.476 


r24n256nw^ 


24 


1.875 


1.5 X lO'^ 


7.6 X 10^ 


4 X 10^ 



















r24n256vzwJ 


24 


1.875 


1.5 X lO'^ 


7.6 X 10^ 


4.7 X 10* 


OC (T 












r24n256cwJ 


24 


1.875 


1.5 X lO'^ 


7.6 X 10^ 


4.7 X 10* 


680 


2 


0.95 


8.7 






r24n256swJ 


24 


1.875 


1.5 X lO'^ 


7.6 X 10^ 


4.7 X 10* 


340 


2 


0.24 


4.4 






r48n256nw 


48 


3.75 


1.2 X 10* 


6.1 X 10* 


~ 3 X IQio 














0.182 





r48n256vzw 


48 


3.75 


1.2 X 10* 


6.1 X 10* 


3.7 X 10^ 


OC (T 


OC (7^^ 


0.52* 


5.1 


0.081 


0.257 


r48n256cw 


48 


3.75 


1.2 X 10* 


6.1 X 10* 


3.7 X 10^ 


680 


2 


0.95 


8.7 


0.046 


0.179 


r48n256sw 


48 


3.75 


1.2 X 10* 


6.1 X 10* 


3.7 X 10^ 


340 


2 


0.24 


4.4 


0.116 


0.448 



Box size in comoving /i~^Mpc. 
* Gravitational softening length in comoving /i~^kpc (equivalent Plummer). 
Galaxy mass resolution limit in Mq; see discussion in il3.1l of no- wind case. 
<* r, = M„i„d/SFR. 
Total feedback energy divided by SN energy (Cliabrier IMF). 
^ Total feedback momentum divided by momentum output by 10^ yr old stellar population (Chabrier IMF), 
f Global mass fraction of baryons in stars at 2 = 0. 
Global mass fraction of baryons that have been ejected in winds by 2 = 0. 
' Volume-averaged to 2 = 0; in this model it is oc a. 
^ Owing to a hardware failure, we only have information at z=0,l,3 for these runs. 



actly the same escape fraction of outflowing material from 
galaxies of all masses. In actuality, hydrodynamic sl owing 
plays an important role l|Qppenheimer fc Dav e 2008), and 
because more massive galaxies live in denser gaseous en- 
vironments, this retards winds more around more massive 
galaxies, causing greater return o f wind material to these 
systems jOppenheimer et "alll2010l ). 

The constant-»7 models (sw and cw) behave more like 
the canonical model for winds, in the sense that at low 
galaxy masses their (constant- velocity) winds can escape the 
potential well, while at high masses they cannot. One can 
envision this as an "effective" mass loading factor that is 
roughly the input va lue of 2 at small masse s, but drops to 
zero at large masses (|Finlator fc Davil2008l ). Alternatively, 
one can view this as rapid re cycling of wind material b ack 
into galaxies at high masses (jOppenheimer et al.ll2010l l. In 
any case, the constant wind speed imparts a feature at a 
mass scale corresponding to where the assumed wind veloc- 
ity is comparable to the galaxy escape velocity. Because of 
the different wind speeds in cw and sw, this feature should 
occur at different masses (roughly a factor of eight apart in 
halo mass, given the factor of two difference in wind speeds). 
Hence comparing these two models enables us to identify fea- 
tures in scaling relations that are related to how the effective 
mass loa ding factor scales with mas s; this is discussed exten- 
sively in lOppenheimer et al. see e.g. their Figure 4). 
As we will argue later, it is the scaling of the (effective) mass 
loading factor with M« that primarily governs most galaxy 
scaling relations. 

Our twelve runs are listed in Tabled They were run at 
the University of Arizona's SGI cluster, ICE, the National 
Center for Supercomputing Applications' Intel cluster, Abe, 
and the University of Massachusetts' Xeon cluster. Eagle. 
Each run took between 100,000 and 150,000 CPU hours, al- 



though nearly 1 million CPU hours were used when includ- 
ing smaller box sizes and various test runs. We note that 
due to technical difficulties, we only have galaxy catalogs 
for the r24 runs at z = 0, 1, 3, and do not have any informa- 
tion regarding satellite galaxies in these runs. However, all 
the important trends can be seen from the other runs at the 
redshifts where the r24 runs are not available. 

We emphasize that all these models are heuristic and 
intended only to explore how these wind scalings impact 
galaxy properties. We are not directly advocating any par- 
ticular physical process, though we are attempting to use 
physically-motivated recipes. In particular, the momentum- 
driven wind scalings can arise from models that only em- 
ploy supernova energy to drive winds, owing to the in- 
teractions of_su£er4Ma^ed_JS^^ with surrounding 
gas (e.g. iDalla Vecchia fc Schavd |2008| ). Hence the name 
"momentum-driven" is only intended to indicate the scal- 
ings rather than a detailed physical model. 

Nevertheless, owing to its prominence in this paper 
and others, we briefly discuss some issues related to the 
physics of the momentum-driven scalings wind model. The 
basic physical motivation is a model in which radiation 



generate an outflow ( 


Murray. Quataert, fc Thompsonll2005l: 


IZhane & Thompson! 


2010l). First, it is important to note 



that the momentum input required to drive winds at hun- 
dreds of km/s with mass loading factors as assumed above 
significantly exceeds the photon momentum emitted by the 
stellar populations (M. Haas fc J. Schaye, priv. comm. iHaad 
2010). We show this ratio in Table[T] an d is a factor o f 5.2 for 
our r48n384vzw simulatiorQ Following iHaad l|2010D . this is 



^ The momentum requirements are similarly large for the 
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calculated from the luminosity of stars for a solar-metallicity 
10^ yr old stellar population assuming a IChabrieil (|2003l ) 
IMF, divided by c. The large value indicates that momen- 
tum alone, under the assumption of single photon scattering, 
cannot power the outflows in the vzw model. 

It is possible that subsequent re-radiation of photons 
could cause significantly more momentum deposition than 
what is emitted by the stars if the gas is optically thick 
to infrared radiation. For an infrared optical depth of tir, 
momentum deposition is increased by 1 + r/_R. This factor 
can only be significantly gre ater than unity in highly dense 
:as w here young stars form. iHopkins. Quataert. fc Murravl 
201l[ l computed rm in simulations of radiation-driven 
winds from high-resolution disk galaxies, and indeed found 
that in some cases tir can approach an order of magni- 
tude, driving an outflow velocity of ~ 200 km/s leaving the 
disk with mass loading factors of order unity. While this 
is encouraging, the issue of whether enough momentum is 
available from stars, and whether it can effectively couple to 
the dust and in turn to the gas, remains an open question. 

Another possibility is that supernovae and photons from 
young stars work together to drive winds (e.g. iNath fc Silkl 
I2OO9I : iMurrav. Menard, fc ThompsonI I2OI0I ) . Invariably, this 
must be happening at some level, s ince both processes 
are energetically comparable. iHaad (|2010l . Figure 2.19) 
showed that the momentum from SNe can be an order 
of magnitude larger than that from stellar radiation at 
10^ yr after the burst, although it dro ps rapidly with time. 
iMurrav. Menard, fc ThompsonI (|2010l ') showed that for a 
300 km/s wind, radiation pressure can yield a mass load- 
ing factor up to 77 = 0.4 in the single-scattering regime, 
which is factors of many less than what our momentum- 
driven wind scalings require. Hence supernova momentum 
input would have to dominate overall, particularly at large 
radii where the optical depth to absorb photons becomes 
small. Despite being driven primarily by supernovae and 
not radiation pressure, the scalings predicted in this sce- 
nario continu e to follow those expected for mom entum- 
driven winds (|Murrav. Menard, fc Thompsonll201Cf ). 

In short, the momentum-driven wind scalings that we 
have assumed challenge the available momentum budget, 
but could plausibly arise from a combination of physical 
processes. Much work remains to be done to marry large- 
scale constraints from cosmological simulations with de- 
tailed studies of wind expulsion from star-forming regions. 



2.3 Galaxy and halo identification 

We use SKIE0 (Spline Kernel Interpolative Denmax) to 
identif y galaxies as boun d groups of star- for ming gas and 
stars (|Keres et al.l l2005l : lOppcnhcim er et al.l [2OIO). Since 
we specifically include just star-forming gas, our simulated 
galaxies only account for gas within the galactic ISM and not 
from the extended halo; this defines the galaxy's gas mass. 
Our galaxy stellar mass limit is set to be 64 star parti- 
cles (Finlator et al. 2006), resulting in the masses listed un- 
der mgai in Table [T] We will argue that for the no- wind case 



constant-Uu, models, but in this case the physical motivation is 
supernova energy, so the momentum input is not relevant. 
^ http://www-hpcc.astro.washington.edu/tools/skid.html 



the resolution convergence is poorer, and we list approxi- 
mate masses that show convergence. Unless otherwise noted, 
we will only consider galaxies with stellar mass M* mgai 
in our analysis. 

We separate galaxies into central and satellite galaxies 
by associating each galaxy with a h alo. We identify h alos 
via a spherical overdensity algorithm (|Keres et al.ll2005l ). In 
brief, we begin at the center of each galaxy and expand out 
spherically until t he mean density enclose d equals that of 
a virialized halo (|Kitavama fc Suto I Il996l ) . We then sub- 
sume all galaxies whose centers are within that of a larger 
galaxy's halo into that larger halo. In this way, each galaxy 
is a member of one and only one halo. The central galaxy 
is the largest one that is at the center of the halo, and any 
other galaxies within that halo are satellites. 



3 STELLAR MASS AND STAR FORMATION 
RATE FUNCTIONS 

3.1 Stellar Mass Function 

We begin by considering the stellar mass function and its 
evolution in our simulations. Stellar mass appears to be 
a governing parameter for galaxy properties, having fairly 
tight scaling relations with many other fundamental galaxy 
properties, re latively independently of environement (e.g. 
iLi et al.llioiol '). Hence the stellar mass function provides a 
key barometer for simulated galaxy populations. 

Figure[T]shows the galaxy stellar mass function (GSMF) 
out to z = 3 in our four ou tfiow simu l ations . Data 
is shown f or comparison from iBeU et all ||2003| . z ~ 
0, poin ts). iBaldrv. Glazebrook. fc Driver! (|2008l . z ~ 0, 
curve), Kaiisawa et al.l ( 20091 . z = l72, 3, curves), and 
iMarchesini et all ( 20091 . z = 2,3, curves). These were se- 
lected as fairly recent determinations of the stellar mass 
function from some of the deepest available near-infrared 
observations, though it is far from a complete sample of 
available data. The curves represent Schechter function best 
fits, while the data points (from different samples) give an 
indica tion of the uncertainties. Note that the . Kaiisawa et al.l 
(|2009h data plotted at 2 = 1 is actually centered at z ~ 0.75, 
but the evolution is not rapid dur i ng th is period. Also, we 
have scaled all data to a lChabrieil l|2003D IMF (as assumed 
in the simulations), by dividing the stellar masses by 1.25 for 
those quoted using a diet Salpeter IMF (n amely, iBell et al] 
l2003l : iBaldrv. Glazebrook. fc DriveJ I2OO8I). and bv 1.8 fo r 



those quo ted using a Salpeter I MF |K ai isawa et al.ll2009l '). 
[Marchcsi ni et al.l |2009i ') used a iKroupa (2001) IMF which 
is very simil ar to Chab r ier, so we apply no IMF correction. 
Finally, the iBell et all l|2003l ) data has been corrected to 
h — 0.7; all other data shown assumed this value for the 
Hubble constant. 

A test of numerical resolution convergence is provided 
by comparing the solid, dotted, and dashed lines for each 
model, which are from the r48n384, r24n256, and r48n256 
scries of runs, respectively. The difference in mass resolu- 
tion between r48n384 and r48n256 is a factor of 2.4, while 
between r48n384 and r24n256 it is 3.3; without accounting 
for this, convergence will appear slightly better at low reso- 
lution. In actuality, for the wind models at every redshift the 
GSMFs show very good convergence (modulo some stochas- 
tic fluctuations at the massive end), and the faint end of 
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Figure 1. Stellar mass functions at 2: = 0, 1, 2, 3 in cosmological hydrodynamic simulations employing our four galactic outflow scalings: 
momentum-conserved winds (green), constant winds (red), no winds (blue), and slow winds (magenta). Solid lines show results from the 
r48n384 series of runs, dotted lines show r24n256 series, and dashed lines show r48n256 series; these illustrate resolution convergence, 
although when interpreting convergence one should keep in mind that these runs are not equally (logarithmically) spaced in mass resolu- 
tion: r48n256 is 2.4x r48n384, which is 3.3x r24n256. The momentum-conserved wind case shows Poisson error bars as a representative 
example. Cyan dashed curve at 2: = 3 shows the result of adding a Gaussian random scatter of cr = 0.3 de x to each galaxy's stellar mass 
in the momentum-conserved win d case. Data po i nts at z = from I Bell et al.l l|2003l l. cur ve at 2: = fromlBaldrv. Glazebrook. fc Driver! 
1I2OO8I I, curves at 2 = 1, 2, 3 from lKaiisawa et al. 1 1I2OO9I '). and data points at 2 = 2, 3 from lMarchesini et al. I l l200g|'l . 



the higher-resolution run joins smoothly onto that of the 
larger-volume run. For the no-wind simulations, this is not 
the case: The higher-resolution run deviates from the lower- 
resolution run at a A/* that is much larger than the nominal 
64 star particle galaxy mass resolution. The reason can be 
traced back to resolution convergence for halos, as opposed 
to galaxies: Halos also require approximately 64 particles 
to robustly host galaxies, and in the no-wind case this cor- 
responds to a much higher M* than in the wind models 
where st ar formation at a given halo mass is much sup- 
pressed (|Oppenheimer et al.ll2010l ). Hence convergence for 
the no- wind case differs from the wind runs, and its resolved 
galaxy mass is higher by roughly an order of magnitude (as 



indicated in Table [T]). The low- mass peak of the GSMF is 
shifted over by roughly the difference in mass resolution be- 
tween the two models, further indicating that the GSMF 
here is impacted by resolution effects. Hence the turn-down 
in the no-wind GSMF at low masses is a numerical arti- 
fact, not a prediction. For the wind models, this turn-down 
happens below the nominal stellar mass resolution. 

The no- wind case (nw; blue), where resolution- 
converged, shows the classic behavior that the GSMF well 
exceeds observations at all masses (see e.g. iKeres et al.l 
l2009bl . for a more detailed discus sion). This is symp- 
tomatic of the ove r cooling problem (|White fc FrenkI Il99ll : 
iBalogh et al] I2OOII : iDave et al.) I2OOII ). in which there are 
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too many baryons locked in stars at all epochs without 
strong feedback. We note that this occurs at all redshifts 
and masses back to z ~ 4 (and ev en farther back; see 
iDave. Finlator. fc Oppenheimer II2006I ). indicating that the 
overcooling problem is not just a late-time phenomenon and 
must be mitigated in the early Universe. 

It has long been suggested that galactic outflows could 
suppress star formation to solve the overcooling problem. 
The magnitude of suppression, factors of several even at 
maximum galaxy formation efficiency around M* increas- 
ing to much larger factors at larger and smaller masses, in- 
dicates that outflows must eject a substantial amount of 
gas from galaxies and/or suppress accretion by adding en- 
ergy to the surrounding IGM. In star-forming galaxies (typ- 
ically Mt < M*), the dominant energy source is star for- 
mation, either from supernovae or young stars. This moti- 
vates all of our wind models having mass loading factors 
exceeding unity. For concreteness. Table [T] shows as its fi- 
nal two columns the global mass fraction of baryons in stars 
at z = and ejected in winds by z = 0, respectively. For 
the momentum-conserved run, the ratio between the two is 
roughly 3.5, i.e. 3.5 x more baryons have been ejected from 
galaxies than have formed into stars by z = 0. The constant- 
?7 models have a ratio of approximately 3.8. That this ratio 
signficantly exceeds unity is a generic feature of all models 
that solve overcooling via an ejective feedback mechanism 
such as galactic outflows. 

The exact shape of the z — GSMF depends sensi- 
tively on outflows, albeit in a subtle way. As discussed in 
lOppenheimer et al.l (|2010l ). it is not the overall suppression 
by ejection that governs the shape, but rather the effective- 
ness of subsequent re-accretion of ejecte d material, i.e. recy- 
cled w ind mode accretion. As shown in lOppenheimer et al.l 
l|20ldl , without wind recycling, all GSMFs look fairly sim- 
ilar, and show a steep power law behavior with a faint-end 
slope comparable to that of the halo mass function. Because 
wind material re-accretes into larger galaxies faster ( "differ- 
ential recycling"), it boosts high-mass galaxy masses rela- 
tive to low-mass ones, and therefore flattens the faint end of 
the stellar mass function. At sufficiently small masses, the 
timescale for wind re-accretion exceeds the Hubble time, at 
which point the mass function returns towards tracking the 
dark matter halo mass function as there is no differential ef- 
fect from recycling. At sufficiently large masses, recycling is 
so fast that it effectively is like having no winds at all. This 
creates a three-tier behavior for the GSMF in simulations: It 
follows the steep dark matter mass function both above and 
below the mass range where differential recycling operates, 
and it is flatter within that mass range. 

All our wind models display this three-tiered behavior, 
though it is less evident in some cases. At high masses, they 
have a slope similar to that of the no- wind model. In the in- 
termediate mass regime, the slope is flatter or even inverted, 
and finally at low masses once again the GSMF is reverts to 
being steep. We note that all models grossly overpredict the 
meissive end of the GSMF compared to observations. This is 
because we do not have any feedback mechanism to truncate 
star formation in massive gala xies, such as feedback f rom ac- 
tive galactic nuclei (AGN; e.g. lDi Matteo et al.ll2005l '). Hence 
our first key point is that our current simulations can only 
reliably probe the regime of star-formi ng galaxies, whic h 
dominate at masses < M* f» lO^M© dSaUmeTal] |20o3) . 



In lGabor etHI (|2010l ) we show that it is possible to include 
empirical feedback mechanisms to quench star formation in 
massive galaxies, without significantly affecting the popula- 
tion of galaxies below Af * . This is consistent with the obser- 
vation that strong AGN activit y is not seen in lower-m ass 
galaxies today (M, < 10^° '^Mm: iKauffmann et~aLll2004 ). 

The exact transitions between the various GSMF tiers 
depend strongly on wind model. In the constant-77 cases, the 
three-tiered behavior is manifested as a bump in the GSMF 
at high masses at z = 0. In the constant wind case (cw; red), 
the middle tier covers a small mass range, and differential 
recycling is so strong that it actually produces an inverted 
slope of the GSMF within that small mass range. Outside 
that, the GSMF is quite steep, roughly following the slope 
of the no-wind case. The slow-wind case (sw; magenta) fol- 
lows a similar pattern, but the mass scale of the differen- 
tial recycling is lower because wind speeds are slower, and 
hence wind recycling becomes effective at a smaller scale. 
In neither case does the GSMF resemble the observed one 
at 2: = even in the regime where quenching feedback is 
unimportant, although the slow-wind case is not as far off. 

The momentum-driven scalings case (vzw; green) pro- 
duces a more gr adual differential recyclin g curve than the 
constant-77 cases (jOppenheimer et al.H2010l ). Hence the mid- 
dle tier is not inverted, but merely shallower than the high 
and low mass ends of the GSMF. This mimics the behavior 
exhibited by the observed GSMF. R ecent observations by 
iBaldrv. Glazebrook. fc Driveij l|2008l ) have conclusively de- 
tected an upturn to a steep faint-end slope. In detail, the 
vzw simulation is still mildly steeper than observed and the 
upturn occurs at slightly too high a mass. If our physical 
interpretation for the origin of the GSMF is correct, this 
would imply that the momentum-conserved wind model has 
a differential recycling curve that is close to, but not quite, 
correct. It is worth mentioning that modeling the dynamics 
of wind recycling is subject to significant num erical difficul- 
ties as discussed in lOppenheimer et all l|2010l , though more 
so at in higher-mass halos where hot gas is present), so the 
level of discrepancy with data is probably within modeling 
uncertainties, but more careful simulations are needed to 
determine this. 

The evolution of the GSMF can now be probed to high 
redshifts thanks to advancing deep near-IR surveys. At the 
massive end, the observed GSMF has a less pronounced cut- 
off at high-z. Since this truncation is associated with the 
presence of "red and dead" galaxies, the relative dearth of 
such galaxies at z > 2 explains why the truncation is much 
less sharp. At z ~ 2 — 3, the vzw run yields a GSMF at 
Af, > W^'^'^ Mq that is in fair agreement with data, al- 
though it appears somewhat too steep. Meanwhile, sw over- 
produces the number of galaxies at these stellar masses 
similar to the no-wind case. The constant wind case gen- 
erally matches well here, although it still shows a bump 
at Mt > 10^^ Mq that is not seen in the data. Consider- 
ing all redshift from z — — 3 in toto, the momentum- 
driven scalings model yields the best match from amongst 
these models to the population of star-forming galaxies at 
10^" < M. < 10^^ Mq. 

The faint-end slope of the GSMF has been observed to 
evolve towards being steeper at high redshifts. Our simu- 
lations follow this trend qualitatively, but only within the 
intermediate mass regime where differential recycling flat- 
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tens (or inverts) the GSMF slope at later times. Outside of 
this regime, the slope always approaches the faint-end halo 
mass slope of a « 2. For instance, in the vzw run, the inter- 
mediate mass regime occurs at 10^'^ < M» < lO^^'^M©, and 
the slope at z = is a ~ 1.4, which is close to the observed 
slope in that mass range. In the sw case, the intermediate 
mass regime is narrower, 10^° < M* < 10^°'^ Mq, while in 
the cw case it is equally narrow but occurs at a mass ~ 8x 
higher. In both constant-77 simulations, the GSMF slope at 
= is inverted in this regime, i.e. a < 1, but at higher 
redshifts a increases. 

The general reason for the increase in a within the 
intermediate mass regime is that wind recycling becomes 
less important at higher redshifts. This is primarily because 
the re is less time for ejected material to return to galax- 
ies l|Qppenheimer et al.ll201ol ). Our simulations, particularly 
the momentum-conserved case, tend to find median recy- 
cling times that are fairly c onstant at ~ 1 — 2 Gyr, in - 
dependent of cosmic epoch (|Oppenheimer fc Pavel |2003 ): 
this is partly because wind speed s at a given M, are 
higher ijOppenheimer fc Pavel I2OO8I ) at high-z. Hence the 
features associated with recycling become less prominent, 
and in fact they are essentially invisible at 2 > 2 in the vzw 
case. 

Obs ervationally, a is seen to evolve as a ~ ao + (0.082 ± 
0.033)z (jFontana et al.ll2006l ). where ao 1.2 - 1.3. Ef- 
fectively, this corresponds to the slope in the mass regime 
of 10^" < A/, < lO"Af0, which is what can be probed 
in practice out to high-z at this time; at lower masses, 
the present-day GSMF becomes steeper, but this cannot 
be reliably traced at high-z. The momentum-driven scal- 
ings model matches well at low- 2; (below M*), but at high-z 
the faint end is quite a bit steeper: e.g. at z = 2, a ~ 2 
in the regime where it is observed to be more like 1.5 or 
so. Hence while the evolutionary trend is qualitatively as 
observed, there are quantitative discrepancies in the sense 
that this model has too many small galaxies (as is the case 
with all models). This has also been noted at 2 ~ 4 — 7 in 
a compariso n of various simul ations to mass function obser- 
vations (Gonzalez et al.ll201ll ): all models have too steep a 
mass function at A/* < Mq. This may indicate that feed- 
back processes in very small systems are not being correctly 
represented in these models, a point that will be reiterated 
throughout this paper. 

One issue that may explain part of the difference be- 
tween the simulations and data are observational uncertain- 
ties in determinations of stellar masses, which can become 
significant particularly at high-z. To explore this, at z = 3 
we added a Gaussian random scatter with u = 0.3 dex to 
each galaxy's stellar mass in the momentum-conserved wind 
simulation (which is a typical observational uncertainty; see 
iMarchesini et all 120091 *) . and recomputed the stellar mass 
function. The result is shown as the dashed cyan line in Fig- 
ure [T] Because the mass function is steep, more small galax- 
ies get scattered to larger masses than vice versa, and the 
mass function does becomes slightly shallower (and in this 
case, agrees slightly better with observations at the massive 
end). But this does not go far towards reducing the overpro- 
duction of small galaxies in the simulations. There are, addi- 
tionally, substantial systematic uncertainti es in determining 
stella r masses from photometric data (e.g. IMarchesini et all 
I2OO9I ): we have not considered such effects here, but it's con- 



ceivable that they could dominate over the statistical errors. 
Hence while this effect should be accounted for when com- 
paring carefully to data, it does not qualitatively alter the 
faint end discrepancy. 

Overall, the momentum-driven scalings model does the 
best job of the models considered in reproducing the ob- 
served GSMFs around Al* from z — 3 0. There are 
still significant discrepancies at higher and lower masses, 
so the agreement is only good in the range of 10^" < Af, < 
10^^ Af0. Nevertheless, this is a significant success, a first for 
cosmological hydrodynamic simulations (that we are aware 
of). While the discrepancy at high masses is well-known and 
strongly suggestive of quenching feedback absent in these 
simulations, the discrepancy at small masses is less cer- 
tain owing to still substantial observational uncertainties at 
high-z. Upcoming deeper observations should quantify this 
discrepancy to greater precision, possibly providing insights 
into the nature of feedback processes in high-redshift dwarf 
galaxies. 

3.2 Star Formation Rate Function 

Figure [2] shows the star formation rate function (<E>sfr) for 
galaxies in our four wind models at four redshifts. The dot- 
ted line shows the SFR at the galaxy stellar mass limit listed 
in Table [Jl obtained by fitting the relationship between SFR 
and M, as displayed in Figure [3| Since we do not have a clean 
stellar mass threshold in this plot, we show all galaxies, but 
$SFR at SFRs below the dotted line may not be robust. Be- 
cause the relationship between SFR and Af, varies slightly 
with wind model, each model has a different SFR threshold. 
Note that the no-wind case has a significantly higher re- 
solved galaxy mass threshold in Table[T]for reasons described 
in the previous section, and hence its dotted line is at a much 
higher mass. O bservations at z = from a GALEX+IRAS 
compilation by iMartin et al.l (|2005l) are shown as the solid 
black line, at z ^ 1 fromlLv et al] 1I2OIOI), and at z = 2 from 
an Ha survey bv lHaves. Schaerer. fc 7)stlinl l|2010l ). 

At z = 0, t&sFR shows different behaviors amongst the 
various wind models that mimic many features seen in the 
GSMF. The effect of differential recycling is seen in the 
three-tiered behavior of the constant-?; models, but the inter- 
mediate regime is not obviously seen in the vzw case. Never- 
theless, relative to no-winds, vzw shifts star formation from 
smaller to larger systems, owing both to its mass loading de- 
pendence as well as differential recyclin g. The large amount 
of recycled wind accretion at late times (jOppenheimer et al.l 
I2OIO) produces a "^sfr at high masses that is actually higher 
for all the wind models than for the no- wind case. The winds 
are metal-enriched and so are particularly effective at cool- 
ing back onto the galaxy (often failing to shock heat to the 
temperature of surrounding halo gas despite their large out- 
ffow velocities), and they delay star formation towards later 
epochs making the effect at high masses stronger towards 
lower redshifts. 

Moving to higher redshifts, the strength of the differen- 
tial recycling features in the cw and sw simulations lessen 
as with the GSMF, but they are still prominent even at 
the highest redshifts. For vzw, "1?sfr becomes quite a bit 
steeper at high-z. The shape of the no-wind $sfr does not 
vary, but the amplitude evolution of $sfr is fairly constant 
from z = 3 2, and then drops with time reflecting the 
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Figure 2. Star formation rate functions •I'sfr at 2 = 0, 1, 2, 3 in our r48n384 series of simulations employing our four galactic outflow 
scalings: momentum-conserved winds (green), constant winds (red), no winds (blue), and slow winds (magenta). Solid lines show results 
from the r48n384 series of runs, dotted lines show r24n256 series (3.3 X better mass resolution), and dashed lines show r48n256 series 
(2.4 X worse mass resolution). The momentum-conserved case shows Poisson error bars as a representative example. Dotted lines indicate 
the SFR at galaxy stellar mass limit (see Table [ij at each redshift, from a linear fit to the A/*— SFR relation for each model; below 
this SFR, the resul ts are less likely to be robust. The thick line at 2 = shows >1?SFR derived from ultraviolet plus infrared data from 
iMartin et al1ll2 005Y while points at 2 = 1 and thick line at z = 2 show $SFR derived from observed Ha luminosity functions at 2 « 0.81 
bv lLv et alT l l201Q'l and at 2 ;^ 2.2 bv lHaves. Schaerer. fc OstlinI ll2010h . respectively. 



overall drop in cosmic SFR. For example, $sFR(10MQ/yr) 
at 2 ~ 2 — 3 is about --^ 10~^'^, dropping to ~ 10^^ at 
2 = 1 and 10~^ at z — 0. In contrast, the momentum- 
conserved wind model shows a less rapidly evolving $sfr, 
as 3>sFR(lOM0/yr)f» —2.3 at all redshifts, although $sfr 
shows a flatter faint-end slope at low redshifts. 

Compared to observations at 2: = (|Martin et ahlfioosh . 
none of the models fare well, as they all overproduce the 
number of star-forming galaxies at higher SFR's. This is 
not surprising, and can be straightforwardly traced to the 
lack of any feedbac k mechanism to qu ench star formation 
in massive galaxies JCabor et al.ll201(]|) as observed in the 
real Universe (e.g. iKauffmann et al.ll2004 ). Note that the 



number of star-forming galaxies must be suppressed down 
to quite moderately star-forming systems. For instance, in 
the momentum-conserved case, $sfr is overproduced at 
SFR> IMq/ji. Therefore galaxies larger than the Milky 
Way must be increasingly quenched, which coincides with 
the discrepancy in the GSMF above M* . The star forma- 
tion rate function therefore provides a strong constraint on 
models for quenching massive galaxies. 

At higher redshifts, we compare to t he Ha luminos - 
ity function deter mined at z = 0.81 by Lv et al.l (|2010l ) 
and at z « 2.2 bv lHaves. Schaerer. fc OstlinI (120101 ). both 
of which prob e down to ^ IMg /yr. We obtain SFR from 
Lhci using the iKennicuttI l| 1998b ) conversion divided by 1.8 
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to correct from Salpeter to Chabrier IMF, namely SFR= 
4.4 X lO'^^Li/cMo/yr. For ILv et all l|2010l ') we show their 
extinction- and in completeness-corrected 2.5a sample. For 
iHaves. SchaererT^ Ostlin (2010) we show their "combined" 
best-fit Schecht er function w h ich is supplemented by bright- 
end data from iGeach et al.l l|2008l ). and correct for 0.977 
magnitudes of extinction {Av = 1.19 at 6563A) indepen- 
dent of Lho ■ At these epochs, quenched massive galaxies are 
increasingly rare, so the comparison should be more mean- 
ingful than at 2 = 0. 

At both epochs, the no-winds run overproduces the 
number of galaxies with SFR< lOOMo/yr, sw overproduces 
galaxies with SFR> lM0/yr, and cw underproduces galax- 
ies with 1 <SFR< lOOM0/yr. The momentum-driven scal- 
ings model fares reasonably well at both redshifts. At 2 ~ 1 
it still slightly overproduces high-SFR galaxies, which is the 
onset of a trend that becomes much more prominent by 
2 = 0. At 2 ~ 2 the bright end matches well, while the faint 
end is modestly steeper th an observed. Note however that 
iHaves. Schaerer. fc 6stiinl (|2010l ') derive a steeper Ha lumi- 
nosity function from their faint-end data alone (a — 1.72 
instead of q = 1.60 for the fit shown here). Furthermore, 
convolving in the uncertainties in M, determinations will 
tend to make the predicted $sfr shallower, as we showed 
with the GSMF in Figure [T] Hence the discrepancies may 
not be significant. 

Overall, the simulated star formation rate functions 
of the constant-r; models display a three-tiered trend aris- 
ing from differential wind recycling as seen in the GSMF. 
The form and evolution of $sfr provides independent con- 
straints on how outflows govern star formation across cosmic 
time. At low-z, all current models produce too many rapidly 
star-forming galaxies, owing to a lack of quenching feedback 
in massive galaxies. At z > 1, the momentum-driven scalings 
simulation predicts $sfr in good agreement with available 
Ha data, while other wind models do not. 



4 SPECIFIC STAR FORMATION RATE 

The specific star formation rate (sSFR) measures the in- 
tensity of ongoing versus past-averaged star formation, so it 
provides an important barometer for how a given galaxy has 
assembled. It also most directly governs the observed colors 
of galaxies. The main sequence of galaxies reflects a rela- 
tion between sSFR and M, , and yields important insights 
into galaxy growth. In this section we study the specific star 
formation rate in our simulations. We will particularly fo- 
cus on how the slowly-evolving equilibrium between inflows, 
star formation, and outflows governs the behavior of sSFR 
as a function of mass and epoch. 



4.1 Specific Star Formation Rate vs. M, 

Figure [3] shows the sSFR of galaxies at z = in our four 
wind models. A running median within mass bins is shown 
in magenta with la variance; this does not include the non- 
starforming galaxies shown along the bottom of the plot, 
which do not substantially change the median except at the 
lowest masses in the constant-?; models (we will return to 
this in ^5]). The points are color-coded by gas fraction within 



each mass bin, blue meaning high /gas (> 0.5a above me- 
dian), red for low, and green in between. Here we deflne gas 



fraction as /g; 



gaB/(Mg 



-I- M*), with Mgas including 



all gaseous phases in the star-forming ISM. Observations 
of galaxies are shown as the thi ck lines from a co mpila- 
tion of GALEX and SDSS data bv lSalim et akl (120071 ): Solid 
line represents the median s SFR for all galaxies classified as 
having no AGN content in a 'Baldwin . Phillips. &: TerlevichI 
(1981) diagram, while the dashed double power-law repre- 
sents those classified as red (NUV— r > 4). The observed 
median for all galaxies drops sharply above > 1O^^M0 as 
red and dead galaxies dominate in that mass range. 

Simulations genericall y produce a fairly tight relation 
between SFR and M. (e.g. lDave et al.ll200ol : iFinlator et~al] 
I2OO6I ). In all cases, the relatively fiat sSFR(M*) implies 
that star formation rate scales roughly linearly with stel- 
lar mass, and at the high-mass end in all cases the scaling 
is sub-linear. The basic similarity indicates that the sSFR 
is not h ighly sensitive to feedback, as h as been noted pre- 
viously (|Davell200^ : iDutton et al.ll2010l '). The basic reason 
is that, to first order, feedback processes reduce SFR and 
M* in conjunction, thereby keeping sSFR roughly similar. 
All models produce massive galaxies that are forming stars 
vigorously, in contrast with observations that indicate most 
galaxies at high masses are p assive, indicating t he need for 
some quenching mechanism (|Gabor et al. I [20 11'). AU mod- 
els also produce a population of low-mass galaxies that are 
not forming stars; these are predominantly satellite galaxies, 
and will be discussed in ^ 

The basic trends in this plot can be understood in the 
context of a balance b etween accretion, outfiow , and star for- 
mation. As shown in IFinlator fc Pavel (|2008l ). galaxies live 
in a slowly-evolving quasi-equilibrium between these quan- 
tities, such that 



+ SFR, 



(4) 



where the first two terms are the mass infiow and outflow 
rates, respectively. We can rewrite this using the deflnition 
of the mass loading factor = Moutfiow/SFR, so that 



SFR = Minflow/(l + »7)- 



(5) 



This equation, while simple, implicitly makes some non- 
trivial assumptions. First, it assumes that the variations 
in the inflow rate are slow compared to the rate of pro- 
cessing gas into stars, so that galaxies are never strongly 
out of equilibrium (e.g. this is not valid during major 
merger events). Second, it assumes that galaxies do not col- 
lect large reservoirs of gas and then process it, but rather 
process gas as it is made available to the ISM (from in- 
flows and outflows); we call this "supply-regulated" star 
formation, to distinguish it from other scenarios where 
gas aggregates in or around galaxie s before being con- 
sumed rapidly (e.g. Eggen. L vnden- Bell. fc Sand age 1963; 



iMaraston et al.l I2OIC ). Third, there is no explicit depen- 
dence on cosmic epoch or environment, as these are as- 
sumed to be of secondary importance in galaxy growth. 
Fourth, it ne glects addit ional rese rvoirs of gas such as stel - 
lar mass loss (ISchave et al . 2010 : Leitner fc Kravtsov1l2010l ) 
and wind recycling ( Oppenheimer et al. 201(J ) that are par- 
ticularly important at late times, and hence this equation is 
expected to be more valid at ea rlier epochs. We will sh ow in 
this paper and Paper II (see also lFinlator fc Davj200"3 ) that 
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Figure 3. Specific star formation rate (sSFR=SFR/M, ) as a function of stellar mass at 2 = in our r48n384 series of simulations 
employing our four galactic outflow scalings; momentum-conserved winds (upper left), constant winds (upper right), no winds (lower 
left), and slow winds (lower right). Points are color-coded by gas fraction within each stellar mass bin: Galaxies with /gas exceeding 
0.5fT above the median are blue, less than O.bcr below the median are red, and in between are green. Galaxies with sSFR= are shown 
as red points along the bottom. A running me dian is shown in magenta, excluding non-starforming galaxies. Thick lines show observed 
GALEX+SDSS best-fit relations llSalim et al.ii2007.) for non-AGN star-forming galaxies (solid) and blue galaxies with NUV — r < 4 
(dashed, broken power law). 



simulated galaxy properties generally follow those expected 
from the equilibrium relation (eq. |4]), although moreso at 
early epochs, suggesting that at least in these models, the 
above assumptions are broadly valid. 

Within the context of this equilibrium model, let us re- 
turn to examining sSFRs. In the no-wind case, = 0, so 
the trend of sSFR(Mt) reflects the inflow rate into galax- 
ies of a given Af, . Owing to the growth of hot gaseous 
halos th at retard inflow in l arger halos (e.g. iKeres et al.l 
[2005 : D ekel fc Birnboimll200d ). the sSFR is smaller at larger 
masses. This is one manifestation of dow nsizing, namely 
star fo rmation rate downsizing as defined in lFontanot et all 
l|2009l ). This model also produces archaeological downsizing, 
which means that larger galaxies have older stellage popu- 



lations (also called "natural downsizing " iDave et~aL I I2OO6I : 
iNeistein. van den Bosch, fc Dekel II2OO6I ) , as we will show in 
>i5.3l Note that these downsizing trends, while qualitatively 
in agreement with data, do not by themselves p roduce mas- 
sive red and dead galaxies (|Gabor et al.|[201ol ). In the no- 
wind case, the slope of sSFR(M*) agrees well with observa- 
tions, but the amplitude is too low by ~ x2. Hence para- 
doxically, while no winds produce far too much global star 
formation (and stellar mass), they produce too little star 
formation at a given stellar mass today. They also produce 
^ xlO mor e stellar mass fo r a gi ven halo mass (|Keres et all 
l2009al : lOp pcnhcimcr c t al][2O10l ). so that the excess global 
star formation arises from too much star formation per halo, 
particularly at high and low masses. 
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The wind models deviate from the no wind case based 
on rj{M,) and wind recychng. In the cw and sw models, 
since rj is constant, in the simplest scenario both SFR and 
M* are reduced equally. This shifts a given galaxy down in 
M* from the no- wind case by a factor of (1 + 77) = 3, when 
the mass is small enough that winds don't recycle (i.e. the 
effective mass loading factor is identical to the input mass 
loading factor). Since the overall slope of the no- wind case is 
negative, this results in a lower sSFR at a given mass. An ad- 
ditional suppression is provided by preve ntive feedback from 
heat added to surrounding gas by winds (|van de Voort et al.l 
which lowers the accretion rate and hence sSFR at 
a given M*. This is particularly noticeable in the cw case 
which adds su bstantial energy to surround ing gas and the 
even the IGM (|Oppenheimer fc Davell200^ ). At the highest 
masses, all winds are gravitationally and hydrodynamically 
trapped, resulting in an effective mass loading factor of 77 ~ 
as discussed in W2.'2l and so also follows the no-wind case 
trend with a slightly higher amplitude. The reason for the 
amplitude difference is that wind recycling removes mass 
from small systems and adds it back into larger systems 
at later epochs. Ifence wind recycling shifts star formation 
from low to high masses, an effect that is stronger with the 
higher wind speeds owing to the increase in the recycling 
time and increasing preventive feedback from added heat. 
Both these models are in poor agreement with observations, 
as they have too low sSFR by up to an order of magnitude 
at the low-mass end, and the observations do not indicate 
any positive-slope feature in the predicted mass range. 

In the case of momentum-conserved winds, rj varies with 
mass. At the lowest masses where there is no wind recycling, 
?7 oc oc M~^^^ cc M^^^^ (roughly, although in detail 
smaller galaxies tend to have lower stellar fractions, making 
the slope slightly shallower), since our wind model ties a 
directly to the galaxy baryon ic mass Mgai as identified b y 
our group finder (see eq. 6 of lOppenheimer fc Davj[200i ). 
Hence the sSFR shows a slope that exceeds the no-wind 
case by roughly 1/3. This by itself would results in a slope 
of around zero. But differential wind recycling produces an 
additional effect, which in this model is more gradual and 
occurs over a wider range in mass than in the constant-r; 
cases. At intermediate masses, the additional fuel from re- 
cycling causes a positive slope for sSFR(Af,), for the same 
reason a s in the c onstant-7 7 case except the effect i s more 
gradual. iFirmani. 'Avila-Reese. fc Rodriguez-Pueblal (|201(]| ') 
also pointed out that wind recycling tends to increase the 
slope of sSFR(M,). At the largest masses, the slope again 
revert s to the no-wind case a s recycling times are very 
short (jOppenheimer et al.ll2010l ). so winds have little effect. 

The gradual rolling of the slope of sSFR(Af,) is also 
seen in observations, but it less pronounced and really only 
ev ident at masses be low what are plotted here (see Figure 16 
of lSalim et ai] 120071 '). Even at M, < 10^ Mr, the observed 
sSFR(Mt) has a negative slope (|Salim et al.ll2007l l. whereas 
this model predicts a positive slope. We will see in Paper II 
that this is also reflected in the relation between gas fraction 
and M», which has a turnover to low masses that is likewise 
absent in observations. This could signify that wind recy- 
cling should be more effective in small systems, or perhaps 
less material should be ejected from these systems in the first 
place. Just as with the faint end of the high-redshift GSMF, 
this highlights that even our most successful model has dif- 
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Figure 4. Evolution of tlic specific star formation rate at a stel- 
lar mass of M, = 10^° Mq from 2 = 6^-0 in our r48n384 series 
of simulations employing our four gal actic outflow scalin gs. Ob- 
servations are shown at z = 0, 1 f rom lElb az et al ] l|2007l ') and at 
z ^ 2 from [Gonzalez et al. The downward arrow on the 

data point at 2 ~ 2 indicates the ~ X 2 lowering of the UV-derived 
SFR based on Herschel data. Dashed lines show t he typical ob- 
served ler variance abo ut the M«— SFR relation llNoeske et al.l 
I2OO7I : iDaddi et ahllioOTh . The black dotted line shows a power- 
law scaling of (1 -|- z)^''^^ as predicted from cold accretion-driven 
galaxy growth. 



ficulty correctly reproducing observations of dwarf galaxies. 
We will return to this point in ^5.31 

An interesting second-parameter trend is indicated by 
the color of the simulated data points. Galaxies with high 
gas fractions at a given M« lie above the mean sSFR rela- 
tion, and those with low gas fractions lie below. This is not 
surprising since the presence of gas drives star formation 
based on our assumed star formation law, but it provides 
an interesting generic prediction (independent of outfiows) 
of this class of models where star formation is driven by ac- 
cretion. Such a trend would not be expected, for example, if 
mergers drove elevated SFRs at a given M*; in that case, it 
is the gas configuration rather than overall gas content that 
drives elevated SFRs, as gas is much m ore concentrated in 
merg ers owing to tidal dissipation (e.g. iMihos fc HernquistI 
Il996l ). In Paper II we will discuss second-parameter depen- 
dences of the relationships between stellar mass, metallicity, 
and star formation rate in much greater depth, all within 
the context of the aforementioned equilibrium model (eq. |3| 
for galaxy growth. 



4.2 Specific Star Formation Rate Evolution 

Figure [4] shows the evolution from z = 6 — >■ of the 
specific star formation rate at a particular stellar mass of 
M. = lO^^M© in our four wind models. This mass is cho- 
sen because observations can directly measure sSFR at this 
mass at all epochs (with small extrapolations), and it is 
well-resolved by simulations but is expected to be negligi- 
bly affected by quenching feedback. In the simulations, this 
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sSFR is computed as the median sSFR within a stellar mass 
bin of ±0.25 dex around 10^° Mq. 

Most models show a characteristic evolution of the 
sSFR which is essentially a power law in {1 + z). The fun- 
damental physics giving rise to this is gas inflow driven by 
growth of structure. In the cold accretion paradigm that cap- 
tures the basic be havior in hydrodynamic simulations (e.g. 
iKereset al.ll2005l ). the amount of gas inflowing into the 
star forming region is correlated with the gas inflowing at 
the halo virial radius, since cold strea ms efficiently chan - 
nel material to the center of the halo (jPekel et al.ll2009l ). 
Feedback processes that strongly impact the surrounding 
gas to retard infall can caus e significant deviations from 
this llvan de Voort et"aLll2010l ). but let us set that aside for 
now to examine the implications of this simplest scenario. 
The amount of gas entering the virial radius can be esti- 
mated by the total mass accretion rate times the baryon 
fraction. The total mass accr etion rate onto ha los has been 
well measured in simulations (jPekel et al.ll2009l ). and scales 
with redshift as (f + zf-'^^. 

The scaling of (f is indicated by the dotted line, 
arbitrarily normalised to the observed 2 = sSFR. The no- 
wind and slow-wind cases follow this power law essentially 
at all redshifts. The no-wind case is easy to understand, 
since there is no additional physics regulating star forma- 
tion beyond accretion. The slow-wind case also follows this 
because it simply ejects a constant two-thirds of the ac- 
creted material. Hence the SFR reduction translates into 
a similar reduction in M, , so the sSFR evolution is nearly 
identical to the no-wind case. In the sw case, two additional 
processes are at work: Feedback heats surrounding gas to 
retard accretion into the ISM, which is countered by addi- 
tional accretion provided by recycled winds. Note that at 
this stellar mass, abo ut half the z = accretion is recy- 
cled winds foppcnhc imer et al.ll201ol '). and both recycling 
and heating increase to low redshifts. The impact of these 
processes does not appear to be large (barring a remarkable 
cancellation), as the net effect is only a slight increase in 
sSFR in the sw model over the no- wind case. 

The constant wind case shows a substantial depar- 
ture from no- winds, with sSFR lower by ~ x2 — 5, in- 
creasing at late times. As discussed previously, this model 
produces significant heating of surrounding gas, substan- 
tially retarding inflow into th e ISM from the virial ra - 
dius (|Oppenheimer et al] l2010l : Ivan de Voort et al] |2010| ) . 
This effect becomes stronger at low-z as feedback en- 
ergy accumulates in the IGM, as evidenced e.g. by a 
much higher fraction of cosmic gas in the Warm-Ho t 
Intergalactic Medium in this model jPave et al.l l2010al l. 
At M, — 10^°Mp), recycling is rare in this wind 
model (|Oppenheimer et al.l |2010| ). so the net effect is to 
strongly lower the sSFR relative to the no winds case. 

In contrast, the momentum-driven wind scalings model 
increases the sSFR relative to no winds. In this case, the 
relationship to no-winds is not so straightforward to under- 
stand, because the mass loading factor depends on galaxy 
mass and epoch. Overall, the amount of energy input in vzw 
is smaller than in cw. Recyling is dominant at late times at 
Mt ~ 10^" M0, but at early epochs it is not infrequent. The 
net effect is larger suppressions of the sSFR at early epochs, 
while at late times the wind recycling at this mass causes the 
sSFR to be higher relative to a simple inflow model extrap- 



olated from high-z. Hence the vzw model yields a shallower 
evolution of sSFR(z) because it suppresses star formation at 
early epochs and later recycles that material into galaxies. 

Observations show a very rapid rise from z = ^ 2, 
and t hen essentially no evolution out to very high red- 
shift s toaddi et al.ll2007l: IStark et al.|[201ol : iGonzalez et all 
i2OI0l). [Bouche et al.l (|2010l ') showed that this is diflicult to 
reproduce in a simple accretion-driven scenario, and they 
argued for a rather radical modiflcation to the galaxy for- 
mation paradigm in which no stars formed in halos with 
virial masses less than fO^^M© (as well as above lO^^M©). 
This strongly suppresses early star formation until the char- 
acteristic mass of star-forming halos exceeds this threshold 
mass, at which point it reverts to the accretion-driven scal- 
ing. Clearly this hypothesis is incorrect today; there is no 
evident feature in the galaxy population around halo masses 
of ~ fO^^M©, unlike at ~ IO^^Mq where many transitions 
occur. It is also difflcult to accomodate at high-z, as smaller 
halos need to form stars to reionise the Universe and form 
observed high-a galaxy populations (|Munoz fc Loebll201lh . 
Nevertheless, while implausible, this scenario correctly em- 
phasizes that some process must suppress high-redshift star 
formation signiflcantly in order to match current sSFR evo- 
lution observations. 

The observations are themselves uncertain to some ex- 
tent, particularly at high redshifts. For instance, the orig- 
inal z ~ 2 data point (black circle) was taken from dust- 
corrected UV- derived star format ion rates of star-forming 
BzK galaxies (|Daddi et al.l |2007| ). Recent Herschel data 
characterizing the infrared continuum of such sources sug- 
gests that t he extinction correct ions were overestimated by 
up to ~ x2 (|Nordon et al.|[2"oiOl ): this modiflcation is shown 
by the red arrow and circle in Figure ID Prior to this correc- 
tion, the data exceeded the accretion rate evolution, mak- 
ing it difflcult to reconcile this evol ution within the cold 
accretion paradigm. This prompted IPav^ (|2008l ) to make 
the radical suggestion that the IMF may evolve with red- 
shift, so that the true SFR's were lower by ~ x3. If the 
_ffersc/ie/-derived correction is right, the z ~ 2 data point 
is now consistent with being powered by ongoing acc retion, 
making it viable within the cold accretion paradigm. iPave I 
(|2010al ) pointed out that there still may be a problem with 
respect to the momentum-driven scalings model; this may 
be particularly true if the corrections at these masses (which 
are not directly constraine d by //erscfeej data) are less than 
a factor of two (see also iNordon et al]|2010l ). But this is 
no longer a fundamental difflculty, and this level of discrep- 
ancy could in principle be alle viated by modiflcat ions of the 
feedback prescription. Indeed, ICenel et al.l 120111 ) flnd that 
doing high-resolution disk simulations using a momentum- 
conserved feedback model with somewhat different parame- 
ters results in galaxies that are in very good agreement with 
(corrected) sSFR data at z ~ 2. 

The slope j3 of the main sequence (deflned by SFRoc 
Mf ) is al so observed to evolve downwards fro m being unity 
at high-z jlStark et al.l[2010l:lLabbe et al.l2006[ ). to ^ « 0.9 at 
z ~ 1-2 (IPaddi ct al. '20071: lElbaz et aklboO?!), to ,9 » 0.7 at 
z < 1 (Noeskc Gt_aL 200i) down to z ~ (IBrinchmann et al.l 
I2OO4I ). There is some sensitivity to selection effects, partic- 
ularly in how one chooses a sample of star-forming galaxies 
at the massive end; the preponderance of low-sSFR mas- 
sive galaxies particularly at low-z pulls the slope down de- 
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pending on how many such systems are included in the fit. 
Hence it is not straightforward to compare to our models 
where we do not reproduce large passive galaxies. In gen- 
eral, our momentum-conserved winds case produces slopes 
of /3 ~ 0.9 — 1, roughly independently of redshift, and the 
no-wind case produces a shallower slope. For the constant-r; 
models /3 is ill-defined as they show a strong three-tiered 
behavior. The trend of the slope becoming more shallow 
with time is therefore not obviously evident even in our fa- 
vored model. While this may be a consequence of the lack 
of quenching feedback in mas sive galaxies , we n ote that 
the semi-analytic model of So merville et al.l (j2008l ) includes 
AGN feedback but still produces a slope around unity for 
the present-day main sequence. 

In summary, the evolution of the specific star forma- 
tion rate encodes information about the feedback processes 
that have regulated star formation up to that epoch. Hence 
this is a critical quantity to measure with upcoming surveys. 
The fundamental driver is the gas accretion rate onto halos, 
which evolves as (1 -I- z)'^'^^ , but feedback can modify this 
by either retarding accretion onto the galaxy or adding re- 
cycled wind accretion. Matching observations that suggest 
a significant departure from this scaling at z > 2 requires 
suppressing star formation at early epochs. Of the feedback 
models examined here, only the momentum-conserved scal- 
ings case works in this sense, though it is still well away 
from the unevolving sSFR seen from 2: ~ 4 — 7. Obser- 
vations at these epochs re main fairly uncertain (see e.g. 
ISchaerer fc de Barrod lioiol ') . so we await more robust con- 
straints from upcoming surveys. 



5 SATELLITE VS. CENTRAL GALAXIES 

The equilibrium paradigm for galaxy evolution introduced 
earlier is appropriate for central galaxies, because cold ac- 
cretion channels material to the centers of halos. Satellite 
galaxies, in contrast, must scavenge their fuel from ambi- 
ent halo gas, and therefore can be impacted by a range of 
processes related to their surrounding environment. Satel- 
lite galaxies could thus display a broadly different class of 
properties than central galaxies, for which environment is 
not a major driver except through wind recycling. Here we 
investigate the properties of satellite vs. central galaxies in 
our simulations, contrasting the statistics we have examined 
earlier between the satellite and central galaxy populations. 
We particularly focus on dwarf galaxies where many envi- 
ronmental processes could play a large role. 



5.1 Satellite Fraction 

Figure [5] shows stellar mass functions as in Figure [1] (black 
curves), subdivided into central galaxies (magenta curves) 
and satellite galaxies (cyan curves). Solid curves show z = 
results, while dashed curves show z = 2. Smaller panel below 
each main panel show the fraction of satellite galaxies as a 
function of stellar mass, also at z = (solid) and z — 2 
(dashed). Only the r48n384 series of runs are shown in this 
plot. 

Satellite galaxy mass functions are, for the most part, 
simply scaled-down versions of the central galaxy mass func- 
tions. They display the same features arising from differen- 



tial wind recycling. In general, the properties of satellites ap- 
pear to be established when they are central galaxies in their 
own halo, and are not dramatically affected by falling into a 
larger halo. This is broadly consistent with trends inferred 
from halo occupation distribution models of galaxy cl uster- 
ing within dark matter halos (e.g. IConrov et al]|2006l ). 

For all models at z — 0, central galaxies dominate by 
number at all resolved masses. This may be surprising, as 
it is often assumed that most small galaxies are satellites, 
and that the fraction of satellites increases strongly to small 
masses. The no- wind model does show an increase in satellite 
fraction to lower masses, but only mildly so. In contrast, the 
wind models show a satellite fraction that is approximately 
constant at one-third of all galaxies for M« < M* . For 
M, > M* the satellite fraction decreases, although it is high 
again at very large meisses - these are satellites in galaxy 
groups where the halo occup ation distribution is steeply ris- 

The s e resu lts agree with those 



groups wnere tne nalo occup 
ing (e.g. lBerlind et al.ll200a). 



obtained bv lvan den Bosch et al] (|2008l . see their Figure 
based on the conditional luminosity function constrained by 
the 2dF Galaxy Redshift Survey, and recent SAMs also yield 
similar trends (Guo & W hite 2010. ). At z = 2, there are even 
fewer satellite galaxies at a given stellar mass, typically only 
around 20% of all galaxies for M* < 10^^ Mq. Hence dwarf 
galaxies are not typically satellites, but rather central galax- 
ies in small halos. 

The relatively constant (and low) satellite fraction has 
a significant implication for galaxy formation. It means that 
that whatever physics is responsible for establishing the 
mass function and other such ensemble-averaged observ- 
ables must p redominantly go vern central galaxies, not satel- 
lites (see also lFontanot et al.|[2009il . It is not possible, for in- 
stance, to alter satellite galaxy physics to reconcile the steep 
faint-end GSMF slope in the constant-77 or no-wind simula- 
tions with observations. This is even more true at higher 
redshifts. Hence the evolution of the ensemble population 
of galaxies is governed by central galaxies, and understand- 
ing galaxy evolution is, to first order, tantamount to under- 
standing how central galaxies form and grow. 



5.2 Quenched Satellites 

Figure [5] shows birthrates for galaxies in our four wind mod- 
els, separated into central (blue) and satellite (red) galax- 
ies. The birthrate is defined as the specific star formation 
rate multiplied by the time over which the galaxy has been 
forming stars; the sSFR value s are n oted on the right axis, 
for comparison with Figure [3] iDavel (|2008l ) showed that the 
timescale for star formation in most galaxies is reasonably 
well approximated by the Hubble time less 1 Gyr, since none 
but the largest galaxies (today) form stars in the first Gyr. A 
running median birthrate and lines enclosing 25% and 75% 
of galaxies in a given mass bin are given by the solid and 
dashed lines, for centrals (magenta) and satellites (cyan). 

Birthrates are typically around unity for all models. 
At the massive end all models show a negative slope in 
birthrates with M, , indicative of archaeological downsiz- 
ing. At small masses, this can invert owing to differential 
wind recycling as discussed in the previous section. Over- 
all, however, star-forming galaxies show generically smooth 
and fairly constant SFHs over much of cosmic time to 
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Figure 5. Four larger panels show stellar mass functions at 2 = (solid lines) and z = 2 (dasiicd lines) for our four r48n384 simulations. 
Black lines are the total mass functions (reproduced from Figure [TJ, which are subdivided into into centrals (magenta) versus satellites 
(cyan). Accompanying smaller panels below each large panel show the fraction of satellite galaxies as a function of A/* at 2: = (solid) 
and z = 2 (dashed). 



2 ~ (|Davell200d ). similar to that inferred for the Milky 
Way. 

There are a significant number of galaxies with no ongo- 
ing star formation (i.e. "quenched"), predominantly at low 
masses. The smaller panels show the fraction of quenched 
central (magenta) and satellite (cyan) galaxies as a func- 
tion of Mt. There are a few small central galaxies with zero 
instantaneous SFR's. It is not clear what these are; one pos- 



sibility is that they are actually quenched satellites whose 
orbits carry them outside the virial radius calculated by our 
spherical overdensity halo finder, and therefore are identified 
as centrals. Another possibility is that quenching is partly 
numerical, since this may impact a satellite galaxy's abil- 
ity to cool gas from a hot halo. Since they are very few in 
number, we will leave a detailed investigation of these pos- 
sibilities for the future. 
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Figure 6. Larger panels show the birthrate (i.e. tsFHSFR/_A/», where tsFR = ^Hubble ~ IGyr) for galaxies from our r48n384 runs at 
2 = color-coded by central (blue) and satellite (red) in our four wind models. Galaxies with SFR= ( "quenched" ) are shown along the 
bottom of the plot. Running median (solid curves), including quenched galaxies, are shown for central galaxies (magenta) and satellites 
(cyan). Dashed above and below show running 75% and 25% percentiles, respectively. Small panels below show the fraction of quenched 
galaxies (having birthrate 6 = 0) in a given mass bin for centrals (magenta) and satellites (cyan). 



The satellite quenched fractions become large towards 
small masses. Below a fewxlO'^'M©, the median sSFR drops 
to zero for satellites in all models. Even above this mass, 
satellites show lower typical sSFR's than central galaxies, 
hence environmental processes can affect even large satellite 
galaxies. The different wind models (including no winds) 
show only minor differences in terms of quenched satellite 



fraction, despite significant differences in sSFR. This indi- 
cates that it is not outfiows that are responsible for quench- 
ing satellites, and instead is some process(es) related to en- 
vironment that depends more on the growth of structure 
common to all models. 

Possible environmental processes that quench satellites 
are ram pressure stripping, tidal harrassment, and strangu- 
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lation (i.e. cutting off of the gas supply). ISimha et al.l (|2009l ) 
showed that, in SPH simulations similar to ours, quenching 
of a galaxy occurs on a ~ 1 Gyr timescale after entering 
another galaxy's halo. The process is gradual because sub- 
halos retain their identity for quite some time after entering 
a larger halo, so the satellite does not immediately see the 
full effect of the hot gas in the larger halo. In the cold accre- 
tion paradigm, the gas cooling along filaments tends to fall 
to the center of the halo, and hence galaxies that turn into 
satellites become disconnected from their feeding filaments. 
But strangulation cannot be the entire story, because the 
gas consumption timescales for these systems are of order 
a few Gyr (as we will show in Paper II), whereas quench- 
ing occurs on a timescale of < 1 Gyr. Hence harrassment 
and stripping must be playing a role. Unfortunately, while 
the exact processes that drive satellite quenching are inter- 
esting, t hey are a great num erical challenge to model prop- 
erly (e.g. lAgertz et al~ll2007^ . and the numerical resolution 
and methodology used here is probably insufficient to make 
robust quantitative statements. 

Although there are seen to be a significant number of 
dwarf galaxies with little or no ongoing star formation, S PSS 
data (|Weinmann et al.ll200^ : Ivan den Bosch et al.ir2008l ') ac- 
tually indicate a trend opposite to that predicted here, 
namely a smaller fraction of red satellites (and centrals) at 
smaller masses. One difference may be that these analyses 
defined red based on a color cut between the red sequence 
and blue cloud, which does not straightforwardly correspond 
to being fully quenched or not; the observed galaxies may 
simply have lower specific SFRs. Interestingly, a recent semi- 
analytic model that is more successful than our simulations 
at matching a wider range of data still suffers from the 
same di screpancy, showing an upturn in red fraction at low 
masses l|Kimm et al.l 120091 ). Hence the qualitative discrep- 
ancy may hint at generic difficulties in modeling environ- 
mental effects in satellites within hierarchical models, and 
bears further investigation. 

In summary, the fraction of satellite galaxies in our 
wind models is roughly constant with mass for M* < M* at 
around 20-30%, declining slightly with redshift. Their mass 
functions show trends similar to centrals, likely reffecting the 
fact that these trends are set when satellites were centrals. 
Satellites are increasingly quenched to smaller masses, re- 
flecting environmental processes that gradually disconnect 
these galaxies from their feeding cold streams. High-mass 
quenching must arise from another source that is currently 
not included such as AGN feedback, but low-mass quenching 
can be understood as a result of environmental effects that 
are already included (at least in principle) in these types of 
models. 

5.3 Stellar Ages and Dwarf Galaxy Evolution 

We have seen that below masses of A/. < 10^° M©, the wind 
models show a turn-down in sSFR and birthrate (and in 
Paper II we will see a similar turn-down in gas fraction), 
for both satellites and centrals. This is in contradiction with 
observations that show no such turn-down. 

To investigate this further, we show in Figure [7] the 
mean (mass-weighted) stellar ages of galaxies as a function 
of M,, for satellites (red) and centrals (blue), in our four 
wind models at 2: = 0. This is essentially an archaeologi- 



cal downsizing plot, showing whether more massive galaxies 
have older stellar populations as generally observed. Run- 
ning medians are shown for centrals (magenta) and satellites 
(cyan). 

All models exhibit archaeological downsizing for galax- 
ies at the massive end. But for the wind models, below 
some mass the trend reverses, and smaller galaxies begin 
to be progressively old er. This i s opposite to the trend 
observed by [Pasguali et al.l (|2010l ) in a group catalog de- 
rived from SDSS (see their Figure 4). For the vzw and sw 
models, this happens at M« ~ lO^^M©; above this mass, 
the ages are in reasonable agreement with data. They also 
show that satellites are slightly older than centrals at a 
given mass, wh i ch is a lso qualitatively as observed, although 
iPasauah etal] (|2010l ) find the median mass-weighted age 
difference at M« = 10^°Mq is ~ 0.1 Gyr (becoming smaller 
to higher masses) while our wind models predict a smaller 
difference. The no- wind case shows an archaeological down- 
sizing trend at all masses. 

The mass where the turn-down happens is exactly the 
mass where the specific star formation rate turns over (Fig- 
ure |S} . We argued previously that this owes to wind recy- 
cling, in particular the lack of it in the smallest systems. 
Hence we begin to build a consistent story for the failure of 
wind models in the dwarf galaxy regime: The smallest galax- 
ies in these models begin forming stars early, and owing to 
ejection of mass they lose their gas and cannot reacquire it 
as fast as larger galaxies. This is not because small galaxies 
are satellites: central galaxies show exactly the same trend, 
and most small galaxies are centrals anyway. 

One possibility is that the simulations do not resolve the 
(supposedly) bursty star formation histories of dwarf galax- 
ies, which could give rise to younger ages. It is true that our 
simulations lack the resolution to model internal resonances 
that can drive gas towards the center of galaxies du ring in- 
teractions, fueling a burst (|Mihos fc Hernauistlll996l ). How- 
ever, recent data h as strongly quest ioned whether bursts are 
common in dwarfs. iLee et al1(|2009l ) used the volume-limited 
11 HUGS sample to show that the burst fraction, defined 
as having birthrate> 2.5, is only 6^2% in dwarfs (typically 
< O.IL*), comprising 23lg*% of that population's star for- 
mation. These numbe rs are similar to what is seen for more 
massive galaxies (e.g. iBrinchmann et al]|2004l ) 

In the small panels in Figure [7] we show the burst 
{b > 2.5) fraction as a function of mass, for centrals (ma- 
genta) and satellites (cyan). It is generally quite low, and in 
our favored vzw model it is mostly independent of mass at 
^ 5% for M* < 10^°'^Mq, in both satellites and centrals. 
In the no wind case, the burst fraction rises dramatically at 
small masses, but mostly in central galaxies not satellites. 
The three-tiered shapes of the birthrate curves in the cw 
and sw model result in odd behaviors in the burst fraction, 
and result in very low burst fractions for small systems. The 
favoured momentum-conserved wind case agrees well with 
observations of burst fractions, although the fraction of star 
formation in bursts is ~ 15% which is « la below that ob- 
served. This model's concordance with data suggests that, 
contrary to popular wisdom, small dwarfs are not signifi- 
cantly burstier compared to larger star- forming systems. 

The higher ages, steep high-jz GSMF, and lower sSFR's 
for dwarfs suggest that wind models should suppress early 
star formation in smaller systems even more than is cur- 
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Figure 7. Larger panels show mean mass- weighted age of stars in galaxies in our r48n384 runs as a function of M*, color-coded by 
central (blue) and satellite (red) in our four wind models at z = 0. Running median (solid curves), including quenched galaxies, are 
shown for central galaxies (magenta) and satellites (cyan). The smaller panels show the fraction of satellites identified as "bursting", i.e. 
with birthrate b > 2.5, as a function of mass, for centrals (magenta) and satellites (cyan). 



rently done in the momentum-conserved wind case, and 
produce a larger reservoir of gas for consumption at later 
epochs. This could arise if mass loading factors are larger at 
early epochs than assumed here, and/or wind speeds are 
reduced at late times such that these galaxies can reac- 
quire their ejected gas. This could also arise if the star 
formation law is different in these systems owing to a 
less efficient conversion of atomic gas into molecular (e.g. 



[Robertson fc Kravtsovl l2008t ) - this would tend to delay 
star formation, providing a younger, more vigorously star- 
forming dwarf population today. Investigating this discrep- 
ancy further is likely to reveal new physical processes that 
are important for the evolution of small galaxies at all cos- 
mic epochs. 
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6 SUMMARY AND DISCUSSION 

Cosmological hydrodynamic simulations are now reaching 
a maturity level such that they can be informatively com- 
pared to observations of galaxies across cosmic time. This 
is a new era in this type of modeling, as past such efforts 
have generally produced quite poor agreement with data, 
and have been severely limited by dynamic range particu- 
larly when evolved to low redshifts. With improvement in 
both computing power and input physics, simulations can 
now plausibly reproduce a wider suite of observed relations 
(though still far from all), which in return can yield insights 
into the governing physics. A particularly important new 
physical aspect has been the incoporation of galactic out- 
flows, originally motivated to explain IGM enrichment, and 
simultaneously having a wide-ranging impact on galaxies' 
stellar, gas, and metal content. 

In this paper we present a study of how the stellar con- 
tent and star formation rates of galaxies are impacted by 
galactic outflows. We concentrate on understanding the un- 
derlying physical drivers of these properties, and how obser- 
vations can enlighten us on the way in which these drivers 
operate in the real Universe. To this end we present some 
comparisons with observations, showing that including out- 
flows with momentum-conserved scalings provide the best 
overall match to the ensemble of observations considered in 
this work, although notable discrepancies remain. Given the 
broader success of this class of models in matching IGM (and 
other galaxy) properties, and that these scalings are consis- 
tent with those observed for outflows, this adds to the grow- 
ing body of evidence that strong and ubiquitous outflows 
with these scalings are an important piece in understanding 
the overall evolution of galaxies. 

The evolution of the stellar component can be broadly 
understood within the context of a cycle of gas inflow and 
outflow between galaxies and the IGM. This differs some- 
what from the traditional view of galaxy formation in which 
halos and their mergers drive galaxy evolution. The cycle 
of baryons can at its most basic level be analytically rep- 
resented by a slowly-evolving equilibrium between inflow, 
outflow, and star formation (eq.|3)), where the outflows gov- 
ern how much of the inflowing material turns into stars. 
The galaxy mass dependence of the outflow rate thus di- 
rectly governs e.g. the GSMF, the specific star formation 
rate of galaxies, and galaxy ages. A particularly key aspect 
is the idea of differential wind recycling, in which material 
ejected returns to galaxies in a mass-dependent manner. 
This drives noticeable features in observable properties of 
galaxies as a function of mass, providing a way to constrain 
how outflows move material within the cosmic ecosystem 
surrounding galaxies. We reiterate that this scenario does 
not invoke the halo virial radius, environment, or mergers 
as playing a governing role in star-forming galaxy evolu- 
tion. The cold streams feeding star-forming galaxies gener- 
ally take no notice of the virial radius, are broadly unaf- 
fected by environment, and while such streams also carry in 
galaxies that mer ge, the mass grow th rate in such mergers 
is sub-dominant ( Dekel et al.ll2009l ). In contrast, the virial 
radius, environment, and mergers seem to play a pivotal role 
in quenching star formation, and more generally in the evo- 
lution of massive halos hosting quiescent galaxies. 



With that framework in mind, we summarize the key 
conclusions of this paper: 

• Galactic outflows are required to suppress stellar mass 
growth at all epochs. The stellar mass function is quite steep 
and nearly a power law at high redshifts. At lower redshifts it 
develops a more pronounced three-tier behavior, where the 
middle tier reflects the stee pness of differential (i.e. mass- 
dependent) wind recycling l|Oppenheimer et al.ll2O10l ). and 
the higher and lower tiers tend towards the slope of the halo 
mass function. The middle tier (typically just below A/*) 
becomes shallower with time owing primarily to wind re- 
cycling, which returns ejected material faster to more mas- 
sive galaxies. Observations of the z « GSMF also show 
a three-tier behavior, and the momentum-conserved wind 
simulation comes closest to matching its behavior (although 
the three tiers are not obviously evident). 

• The star formation rate functions show many of the 
same general features as the GSMF, but the detailed evo- 
lution, mass dependence, and dependence on outflows are 
somewhat different. Comparisons of these models to present- 
day data are hampered by their lack of quenching feedback 
in massive galaxies in our models. At 2; > 1 where massive 
quenched galaxies are less common, momentum-conserved 
winds produces a good match to the observed star forma- 
tion rate function derived from Ha data. 

• The specific star formation rate's dependence on mass 
is governed by how the effective mass loading factor of 
outflows varies with stellar mass. The no-wind case shows 
too little star formation at a given M*, despite having far 
too much star formation globally. The momentum-conserved 
wind model matches observations well at ~ 0.1 — 1 M* , 
but all wind models show too little star formation at the 
smallest masses. This may indicate that there should be less 
mass loading or more wind recycling at the M* < O.IM* 
compared to what is currently in these models, or else the 
conversion of gas into stars is not being modeled properly 
in low-mass systems. 

• The evolution of the sSFR is generally well-described 
as being driven by cosmic mass accretion into halos, which 
scales as {1 + z)'^'^^ . The no- wind and slow- wind cases follow 
this trend at all redshifts. Observations do so at 2 < 2, but 
at 2 > 2 they depart from this scaling, indicating that small 
high-redshift galaxies must be significantly suppressed. The 
momentum-conserved wind model qualitatively yields this, 
but not enough to match current high-z observations. Pin- 
ning down the sSFR evolution at z > 2 will have a major 
impact on understanding early galaxy evolution and feed- 
back processes. 

• Satellite galaxies are the minority of galaxies at all 
masses probed by these simulations (down to Af, ~ Mq), 
with a fraction of around one-third at z = independent of 
Af« . Hence understanding the growth of galaxies must focus 
on developing a model for central galaxies. The satellites' 
GSMF follows the same trends as that of centrals. How- 
ever, there are many more non-star forming satellites than 
centrals, particularly at lower masses, owing to environmen- 
tal effects wherein the satellites no longer effectively receive 
inflow from the IGM; this trend may be contrary to that ob- 
served. Small galaxies are at most only mildly more bursty 
than larger ones, and the burst fraction is small in accord 
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with observations. Our results challenge conventional no- 
tions that most small galaxies are bursty satellites. 

Taking a broader view, it has long been recognised that 
the gala:xy formation process must provide some feedback 
mechanism that preferentially suppresses mass growth at 
both high and low masses. In this paper we have advocated 
that it is galactic outflows that are primarily responsible for 
regulating the low-mass end. The advantages of invoking this 
particular mechanism are that (i) outflows are observed and 
seem to follow scalings that yield many desirable properties 
in models; (ii) they concurrently enrich the IGM at all cos- 
mic epochs as observed; (iii) they lower the baryon con tent of 
halos preferentially to lower masses as observed (e.g. iDave I 
l2010bl ): and (iv) they are an ejective feedback mechanism 
that works in conjunction with the cold accretion paradigm, 
in which it is difficult to prevent filaments from channeling 
gas efficiently into galaxies. There are other ideas to sup- 
press low-mass star formation by regulating it as a function 
of mass internally within the ISM (i.e. by varying the star 
formation efficiency as a function of M,), or by preventing 
material entering the halo from accreting into small galax- 
ies. But these scenarios must explain IGM enrichment sep- 
arately, and must hide most of the gravitationally-accreted 
baryons in halos in some undetected form which is becom- 
ing increasingly diffi cult to accomodate fe.g. lMcGaugh et al.l 
l20ld : lDai etaUlioid V 

While attractive, the idea of invoking outflows is not 
without its difficulties. For instance, the amount of mass be- 
ing driven out of these galaxies is large, globally many times 
the amount of mass forming into stars. This is a generic 
result, because some process must suppress st ellar masses 
by such large factors at the low mass end (|Keres et al.l 
l2009al :IOppcnhci mer et al.ll2010l ): in our models, this process 
is galactic outfiows, as opposed to ISM physics or preven- 
tive feedback. The physical mechanism(s) that would drive 
such a sizeable amount of material out of galaxies is un- 
clear. The canonical idea of driving outflows from overpres- 
surized bubbles in the ISM tends to (when modeled) cre- 
ate holes by which energy and metals , but little mass, are 
ejected (e.g. iMac Low fc Ferraralll99^ ). and there are diffi- 
culties with entraining cold clouds as observed without dis- 
rupting them. The scalings preferred by our current compar- 
isons to data are those expected for radiation-driven winds, 
but it is not well known how far out the photons' momentum 
can drive dust, and how f ar out the dust remains coupled 
to the gas. For instance, lYoshida. Kawabata. fc Ohivamal 
(|2010l ) use spectropolarimetry to determine that the dust 
in M82 is moving outwards much more slowly than the gas, 
suggesting that the gas and dust are not strongly coupled. 
Furthermore, the momentum budget of typical stellar popu- 
lations is, under the assumption of single photon scattering, 
insufficent to drive outflows as appear to be necessary, in- 
dicating that supernovae or other sources of energy are still 
required. Without a solid understanding of the dynamics 
driving outflows, simulations such as the ones presented here 
cannot be considered a complete description of the physics 
of galaxy formation. Fortunately, observations of outflows 
across cosmic time are gaining substantial t raction (|Martinl 
I2OO5I : IWeiner et al.|[200g : ISteidel et al.ll2O10D . providing em- 
pirical constraints for models while the theory of outflow 
propagation develops. 



Our momentum-conserved winds provides the best 
match to data on star-forming galaxies of the four mod- 
els considered here. But the agreement is only good 
in the rather narrow mass range from 10^" < < 
10"Mo. At lar ger masses, the current models require 
some additional quenching mechanism perhaps associated 



with AGN 



e^ 



Pi Matteo et al.ll2005l : ICroton et al] l2006l : 



iGabor et al. I2OI1I '). At smaller masses, this model produces 



too vigorous early star formation in dwarfs, resulting in too- 
old stellar populations, too steep mass functions, and in- 
sufficient s tar formation today (a result also seen in many 
SAMs, e.g. iFontanot et al.ll2009l : IGuo fc Whitell2010l ). This 
may indicate that feedback in these dwarf systems evolves 
in some way not captured by this particular outflow pre- 
scription, or else that there is some additional physical pro- 
cess that is important. It is possible that the star forma- 
tion law is different in s maller systems to make gas c on- 
sumption less efficient (e.g. iRobertson fc Kravtsovll2008l ). or 
else that these simulations have too low a density thresh- 
old for star formation (owing to resolution limitations) such 
that they do not properly represent star formation in low 
surfa ce brightness, gas-rich systems (e.g. ,Govcrnat o et al.] 
More broadly, none of our simulations can simulta- 
neously reproduce the shallow stellar mass function, archae- 
ological downsizing, and star formation rate downsizing as 
observed; this remains a fundamental challenge for galaxy 
formation models. Even though the concordant mass range 
of the momentum-conserved wind model is fairly narrow, 
it is still a significant step forward as it provides physical 
insights into the processes that govern galaxy evolution at 
all masses, pinpointing both successes and failures. And it 
is worth recognizing that this particular dex in stellar mass 
happens to be where most of cosmic star formation occurs, 
so it is a critical mass range for understanding global stellar 
mass growth. 

Overall, the general simulation-inspired framework for 
how inflows and outflows govern galaxy formation is com- 
pelling from both an observational and theoretical perspec- 
tive, and is well-situated within the current hierarchical 
structure formation paradigm. This makes it an attractive 
scenario to pursue both analytically and numerically. In Pa- 
per II we will examine how this framework provides intuition 
about the gaseous and metal content of galaxies, thereby 
together covering the primary constituents that determine 
the observable multi-wavelength properties of galaxies at all 
cosmic epochs. 
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